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Abstract 

Several km-scale gravitational-wave detectors have been constructed world wide. These 
instruments combine a number of advanced technologies to push the limits of precision length 
measurement. The core devices are laser interferometers of a new kind; developed from the 
classical Michelson topology these interferometers integrate additional optical elements, which 
significantly change the properties of the optical system. Much of the design and analy- 
sis of these laser interferometers can be performed using well-known classical optical tech- 
niques, however, the complex optical layouts provide a new challenge. In this review we give 
a textbook-style introduction to the optical science required for the understanding of modern 
gravitational wave detectors, as well as other high-precision laser interferometers. In addition, 
we provide a number of examples for a freely available interferometer simulation software and 
encourage the reader to use these examples to gain hands-on experience with the discussed 
optical methods. 



1 Introduction 



1.1 The scope and style of the review 

The historical development of laser interferometers for application as gravitational-wave detec- 
tors |47] has involved the combination of relatively simple optical subsystems into more and more 
complex assemblies. The individual elements that compose the interferometers, including mirrors, 
beam splitters, lasers, modulators, various polarising optics, photo detectors and so forth, are 
individually well described by relatively simple, mostly-classical physics. Complexity arises from 
the combination of multiple mirrors, beam splitters etc. into optical cavity systems, have narrow 
resonant features, and the consequent requirement to stabilise relative separations of the various 
components to sub-wavelength accuracy, and indeed in many cases to very small fractions of a 
wavelength. 

Thus, classical physics describes the interferometer techniques and the operation of current 
gravitational-wave detectors. However, we note that at signal frequencies above a couple of hun- 
dreds of Hertz, the sensitivity of current detectors is limited by the photon counting noise at the 
interferometer readout, also called shot-noise. The next generation systems such as Advanced 
LIGO ^MlIS], Advanced Virgo @] and LCGT [37] are expected to operate in a regime where the 
quantum physics of both light and mirror motion couple to each other. Then, a rigorous quantum- 
mechanical description is certainly required. Sensitivity improvements beyond these 'Advanced' 
detectors necessitate the development of non-classical techniques. The present review, in its first 
version, does not consider quantum effects but reserves them for future updates. 

The components employed tend to behave in a linear fashion with respect to the optical field, 
i.e., nonlinear optical effects need hardly be considered. Indeed, almost all aspects of the design 
of laser interferometers are dealt with in the linear regime. Therefore the underlying mathematics 
is relatively simple and many standard techniques are available, including those that naturally 
allow numerical solution by computer models. Such computer models are in fact necessary as the 
exact solutions can become quite complicated even for systems of a few components. In practice, 
workers in the field rarely calculate the behaviour of the optical systems from first principles, but 
instead rely on various well-established numerical modelling techniques. An example of software 
that enables modelling of either time-dependent or frequency-domain behaviour of interferometers 
and their component systems is Finesse [231 [20] ■ This was developed by one of us ( AF) , has 
been validated in a wide range of situations, and was used to prepare the examples included in the 
present review. 

The target readership we have in mind is the student or researcher who desires to get to grips 
with practical issues in the design of interferometers or component parts thereof. For that reason, 
this review consists of sections covering the basic physics and approaches to simulation, intermixed 
with some practical examples. To make this as useful as possible, the examples are intended to 
be realistic with sensible parameters refiecting typical application in gravitational wave detectors. 
The examples, prepared using Finesse, are designed to illustrate the methods typically applied in 
designing gravitational wave detectors. We encourage the reader to obtain Finesse and to follow 
the examples (see Appendix [A| . 

1.2 Overview of the goals of interferometer design 

As set out in very many works, gravitational-wave detectors strive to pick out signals carried by 
passing gravitational waves from a background of self-generated noise. The principles of operation 
are set out at various points in the review, but in essence, the goal has been to prepare many 
photons, stored for as long as practical in the 'arms' of a laser interferometer (traditionally the 
two arms are at right angles), so that tiny phase shifts induced by the gravitational waves form 
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as large as possible a signal, when the light leaving the appropriate 'port' of the interferometer is 
detected and the resulting signal analysed. 

The evolution of gravitational-wave detectors can be seen by following their development from 
prototypes and early observing systems towards the Advanced detectors which are currently in 
the final stages of planning or early stages of construction. Starting from the simplest Michelson 
interferometer [TH] , then by the application of techniques to increase the number of photons stored 
in the arms: delay lines [35], Fabry-Perot arm cavities [T71[TH] and power recycling [12). The final 
step in the development of classical interferometry was the inclusion of signal recycling [321 ISI] > 
which, among other effects, allows the signal from a gravitational-wave signal of approximately- 
known spectrum to be enhanced above the noise. 

Reading out a signal from even the most basic interferometer requires minimising the cou- 
pling of local environmental effects to the detected output. Thus, the relative positions of all 
the components must be stabilised. This is commonly achieved by suspending the mirrors etc. 
as pendulums, often multi-stage pendulums in series, and then applying closed-loop control to 
maintain the desired operating condition. The careful engineering required to provide low-noise 
suspensions with the correct vibration isolation, and also low-noise actuation, is described in many 
works. As the interferometer optics become more complicated, the resonance conditions, i.e., the 
allowed combinations of inter-component path lengths required to allow the photon number in the 
interferometer arms to reach maximum, become more narrowly defined. It is likewise necessary 
to maintain angular alignment of all components, such that beams required to interfere are cor- 
rectly co-aligned. Typically the beams need to be aligned within a small fraction (and sometimes 
a very small fraction) of the far-field diffraction angle, and the requirement can be in the low 
nanoradian range for km-scale detectors [3SJ[25j. Therefore, for each optical component there is 
typically one longitudinal (i.e., along the direction of light propagation), plus two angular degrees 
of freedom (pitch and yaw about the longitudinal axis) . A complex interferometer can consist of 
up to around seven highly sensitive components and so there can be of order 20 degrees of freedom 
to be measured and controlled [31 [57] . 

Although the light fields are linear, the coupling between the position of a mirror and the 
complex amplitude of the detected light field typically shows strongly nonlinear dependence on 
mirror positions due to the sharp resonance features exhibited by cavity systems. However, the 
fields do vary linearly or at least smoothly close to the desired operating point. So, while well- 
understood linear control theory suffices to design the control system needed to maintain the 
optical configuration at its operating point, bringing the system to that operating condition is 
often a separate and more challenging nonlinear problem. In the current version of this work we 
consider only the linear aspects of sensing and control. 

Control systems require actuators, and those employed are typically electrical-force transducers 
that act on the suspended optical components, either directly or - to provide enhanced noise 
rejection - at upper stages of multi-stage suspensions. The transducers are normally coil-magnet 
actuators, with the magnets on the moving part, or, less frequently, electrostatic actuators of 
varying design. The actuators are frequently regarded as part of the mirror suspension subsystem 
and are not discussed in the current work. 

1.3 Overview of the physics of the primary interferometer components 

To give order to our review we consider the main physics describing the operation of the basic 
optical components (mirrors, beam splitters, modulators, etc.) required to construct interferom- 
eters. Although all of the relevant physics is generally well known and not new, we take it as a 
starting point that permits the introduction of notation and conventions. It is also true that the 
interferometry employed for gravitational- wave detection has a different emphasis than other inter- 
ferometer applications. As a consequence, descriptions or examples of a number of crucial optical 
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properties for gravitational wave detectors cannot be found in the literature. The purpose of this 
first version of the review is especially to provide a coherent theoretical framework for describing 
such effects. With the basics established, it can be seen that the interferometer configurations that 
have been employed in gravitational-wave detection may be built up and simulated in a relatively 
straightforward manner. 

As mentioned above, we do not address the newer physics associated with operation at or 
beyond the standard quantum limit. The interested reader can begin to explore this topic from 
the following references. 

• The standard quantum limit [TDJ 

• Squeezing [5^ [55] 

• Quantum nondemolition interferometry [51 125) 

These matters are to be included in a future revision of this review. 



1.4 Plane- wave analysis 

The main optical systems of interferometric gravitational-wave detectors are designed such that 
all system parameters are well known and stable over time. The stability is achieved through a 
mixture of passive isolation systems and active feedback control. In particular, the light sources are 
some of the most stable, low-noise continuous-wave laser systems so that electromagnetic fields can 
be assumed to be essentially monochromatic. Additional frequency components can be modelled as 
small modulations (in amplitude or phase). The laser beams are well coUimated, propagate along 
a well-defined optical axis and remain always very much smaller than the optical elements they 
interact with. Therefore, these beams can be described as paraxial and the well-known paraxial 
approximations can be applied. 

It is useful to first derive a mathematical model based on monochromatic, scalar, plane waves. 
As it turns out, a more detailed model including the polarisation and the shape of the laser beam 
as well as multiple frequency components, can be derived as an extension to the plane- wave model. 
A plane electromagnetic wave is typically described by its electric field component: with Eq as the 



k gives the direction of the wave with k = oj/c 



(^direction of polarisation \ phase offset 
E{x,y,z,t) =(EdK^cos[(ujh -(kf)+(ipj 



Qield amplitude^ a; = 27r / is the angular frequency 



Figure 1: 



(constant) field amplitude in V/m, Bp the unit vector in the direction of polarisation, such as, for 
example, Cy for J>*-polarised light, uj the angular oscillation frequency of the wave, and k = Cfew/c 
the wave vector pointing the in the direction of propagation. The absolute phase ip only becomes 
meaningful when the field is superposed with other light fields. 

In this document we will consider waves propagating along the optical axis given by the 2-axis, 
so that kr — kz. For the moment we will ignore the polarisation and use scalar waves, which can 
be written as 

E{z,t) = Eocos{ujt — kz + (fi). (1) 
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Further, in this document wc use complex notation, i.e.. 



E = ^{E'} 



with 



E' = E', 



iQ exp (i [wt — kzj). 



(2) 



This has the advantage that the scalar amphtude and the phase ip can be given by one, now complex, 
amplitude E'q = £'oexp(i(^). We will use this notation with complex numbers throughout. For 
clarity we will simply use the unprimed letters for the auxiliary field. In particular, we will use 
the letter E and also a and h to denote complex electric-field amplitudes. But remember that, for 
example, in = £^0 exp(— i kz) neither E nor Eq are physical quantities. Only the real part of E 
exists and deserves the name field amplitude. 

1.5 Frequency domain analysis 

In most cases we are either interested in the fields at one particular location, for example, on the 
surface of an optical element, or we want to know the fields at all places in the interferometer but 
at one particular point in time. The latter is usually true for the steady state approach: assuming 
that the interferometer is in a steady state, all solutions must be independent of time so that we 
can perform all computations at t = without loss of generality. In that case, the scalar plane 
wave can be written as 



The frequency domain is of special interest as numerical models of gravitational-wave detectors 
tend to be much faster to compute in the frequency domain than in the time domain. 



E = Eo ex.p{—i kz) . 



(3) 
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2 Optical Components: Coupling of Field Amplitudes 



When an electromagnetic wave interacts with an optical system, all of its parameters can be 
changed as a result. Typically optical components are designed such that, ideally, they only affect 
one of the parameters, i.e., either the amplitude or the polarisation or the shape. Therefore, it is 
convenient to derive separate descriptions concerning each parameter. This section introduces the 
coupling of the complex field amplitude at optical components. Typically, the optical components 
are described in the simplest possible way, as illustrated by the use of abstract schematics such as 
those shown in Figure [2] 
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Figure 2: This set of figures introduces an abstract form of illustration, which will be used in this 
document. The top figure shows a typical example taken from the analysis of an optical system: 
an incident field Ei^ is reflected and transmitted by a semi-transparent mirror; there might be the 
possibility of second incident field Ein2 ■ The lower left figure shows the abstract form we choose 
to represent the same system. The lower right figure depicts how this can be extended to include 
a beam splitter object, which connects two optical axes. 



2.1 Mirrors and spaces: reflection, transmission and propagation 

The core optical systems of current interferometric gravitational interferometers are composed of 
two building blocks: a) resonant optical cavities, such as Fabry-Perot resonators, and b) beam 
splitters, as in a Michelson interferometer. In other words, the laser beam is either propagated 
through a vacuum system or interacts with a partially-refiecting optical surface. 

The term optical surface generally refers to a boundary between two media with possibly 
different indices of refraction n, for example, the boundary between air and glass or between two 
types of glass. A real fused silica mirror in an interferometer features two surfaces, which interact 
with a reflected or transmitted laser beam. However, in some cases, one of these surfaces has been 
treated with an anti-reflection (AR) coating to minimise the effect on the transmitted beam. 

The terms mirror and beam splitter are sometimes used to describe a (theoretical) optical 
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surface in a model. We define real amplitude coefficients for reflection and transmission r and 
t, with < r, i < f, so that the field amplitudes can be written as The 7r/2 phase shift upon 



E2 = rE:i + i tEi 
E4 = rEi+i tE3 



mirror 

(optical surface described by coefficients r and t) 
El E2 

optical axis 



E. 



E. 



Figure 3: 



transmission (here given by the factor i ) refers to a phase convention explained in Section 2.4 

The free propagation of a distance D through a medium with index of refraction n can be 
described with the following set of equations: In the following we use n = 1 for simplicity. 



E2 = El exp(— i k n D) 
E4 = i?3 exp(— i k n D) 



space 

(propagation defined by coefficients D and n) 
El , E2 

*' optical axis 



Ei 



Figure 4: 

Note that we use above relations to demonstrate various mathematical methods for the analysis 
of optical systems. However, refined versions of the coupling equations for optical components, 
including those for spaces and mirrors, are also required, see, for example. Section [2.6| 
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Figure 5: Simplified schematic of a two mirror cavity. The two mirrors are defined by the amplitude 
coefficients for reffection and transmission. Further, the resulting cavity is characterised by its 
length D. Light field amplitudes are shown and identified by a variable name, where necessary to 
permit their mutual coupling to be computed. 



2.2 The two- mirror resonator 

The linear optical resonator, also called a cavity is formed by two partially-transparent mirrors, 
arranged in parallel as shown in Figure [5] This simple setup makes a very good example with which 
to illustrate how a mathematical model of an interferometer can be derived, using the equations 
introduced in Section [2T1 
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The cavity is defined by a propagation length D (in vacuum), the amplitude reflectivities ri, 
r2 and the amplitude transmittances ii, t2. The amplitude at each point in the cavity can be 
computed simply as the superposition of fields. The entire set of equations can be written as 



exp(— ifcZ3) ai 



0-2 

as 
04 



it2a[ 
r2a'i 

exp(— ifcZ3) as 
riflo + iiiag 



(4) 



The circulating field impinging on the first mirror (surface) a'^ can now be computed as 



exp(— ifcZ?) 03 = exp(— i/cZ?) r2a'^ = exp(— i2fcD) r2ai 
ex.p{—i2kD) r2 (i^iao + »'ia3). 



This then yields 



03 = flo: 



ir2ti cxp{—i2kD) 



1 — rir2 exp(— i 2kD) 
We can directly compute the reflected field to be 



04 



r2tlexp{-'i2kD) \_ /ri - r2(rf + ) exp(-i 2fcD) \ 
' " " ~" " 1 — rir2 exp(— i 2fc_D) / 



1 — rir2 exp(— i 2kD) J 
while the transmitted field becomes 

-tit2 cxp(— i kD) 



02 = flo: 



(5) 



(6) 



(7) 



(8) 



1 — rir2 exp(— i 2kD) 

The properties of two mirror cavities will be discussed in more detail in Section |5.1[ 
2.3 Coupling matrices 

Computations that involve sets of linear equations as shown in Section |2.2| can often be done or 
written efficiently with matrices. Two methods of applying matrices to coupling field amplitudes 
are demonstrated below, using again the example of a two mirror cavity. First of all, we can rewrite 
the coupling equations in matrix form. The mirror coupling as given in Figure [3] becomes 

mirror (r, t) 



ai 
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04 



i t r 
r i t 
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Figure 6: 
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and the amplitude coupling at a 'space', as given in Figure |4j can be written as In these 
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04 



exp(— i kD) 
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Figure 7: 

examples the matrix simply transforms the 'known' impinging amplitudes into the 'unknown' 
outgoing amplitudes. 



Coupling matrices for numerical computations An obvious application of the matrices 
introduced above would be to construct a large matrix for an extended optical system appropriate 
for computerisation. A very flexible method is to setup one equation for each field amplitude. The 
set of linear equations for a mirror would expand to 
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(9) 



where the input vectoi[^ Sinput has non-zero values for the impinging fields and Osoi is the 'solution' 
vector, i.e., after solving the system of equations the amplitudes of the impinging as well as those 
of the outgoing fields are stored in that vector. 

As an example we apply this method to the two mirror cavity. The system matrix for the 
optical setup shown in Figure [5] becomes 
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(10) 



This is a sparse matrix. Sparse matrices are an important subclass of linear algebra problems 
and many efficient numerical algorithms for solving sparse matrices are freely available (see, for 
example, [Sj). The advantage of this method of constructing a single matrix for an entire optical 
system is the direct access to all field amplitudes. It also stores each coupling coefficient in one 
or more dedicated matrix elements, so that numerical values for each parameter can be read out 
or changed after the matrix has been constructed and, for example, stored in computer memory. 
The obvious disadvantage is that the size of the matrix quickly grows with the number of optical 
elements (and with the degrees of freedom of the system, see, for example. Section [7|. 

Coupling matrices for a compact system descriptions The following method is probably 
most useful for analytic computations, or for optimisation aspects of a numerical computation. 



^In many implementations of numerical matrix solvers the input vector is also called the right-hand side vector. 
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The idea behind the scheme, which is used for computing the characteristics of dielectric coatings 
[251 HT| and has been demonstrated for analysing gravitational wave detectors [M] , is to rearrange 
equations as in Figure [6] and Figure [7] such that the overall matrix describing a series of components 
can be obtained by multiplication of the component matrices. In order to achieve this, the coupling 
equations have to be re-ordered so that the input vector consists of two field amplitudes at one 
side of the component. For the mirror, this gives a coupling matrix of 



0.4 J t \ -r + t'^ J \ as 



(11) 



In the special case of the lossless mirror this matrix simplifies as we have r^ + i^ = R + T = 1. 
The space component would be described by the following matrix: 



oi \ _ / exp{ikD) \ f 02 

04/1 exp(— ifc£>) y I 03 



(12) 



With these matrices we can very easily compute a matrix for the cavity with two lossless mirrors 
as 



-^cav — Afmirrorl X -Mspacc X A'/jnirror2 (1^) 

— 1 / — rir2e~ — r2e+ + Tie" 
tit2 \ ~'i"2e^ + rie+ - rir2e+ 



(14) 



with e+ = cxp(ifcl?) and e — exp(— i/cZ?). The system of equation describing a cavity shown in 
Equation Q can now be written more compactly as 



oo \ _ -1 / e+ - rir2e -r2e+ + rie \ / 02 
04 / ii<2 V -?'2e^+rie+ - rir2e+ / \ 

This allows direct computation of the amplitude of the transmitted field resulting in 



(15) 



-tit2 cxp{-i kD) 
1 — rir2 cxp(— 1 2kU) 

which is the same as Equation ([8|. 

The advantage of this matrix method is that it allows compact storage of any series of mirrors 
and propagations, and potentially other optical elements, in a single 2x2 matrix. The disadvan- 
tage inherent in this scheme is the lack of information about the field amplitudes inside the group 
of optical elements. 



2.4 Phase relation at a mirror or beam splitter 

The magnitude and phase of refiection at a single optical surface can be derived from Maxwell's 
equations and the electromagnetic boundary conditions at the surface, and in particular the con- 
dition that the field amplitudes tangential to the optical surface must be continuous. The results 
are called FresneVs equations [34]. Thus, for a field impinging on an optical surface under normal 
incidence we can give the reflection coefficient as 

ni-n2 , , 

, , 17) 

rti + n2 

with ni and n2 the indices of refraction of the first and second medium, respectively. The 
transmission coefficient for a lossless surface can be computed as = 1 — r^. We note that 
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anti-reflective 
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optical axis 




Figure 8: This sketch shows a mirror or beam sphtter component with dielectric coatings and 
the photograph shows some typical commercially available examples |13) . Most mirrors and beam 
splitters used in optical experiments are of this type: a substrate made from glass, quartz or fused 
silica is coated on both sides. The reflective coating defines the overall reflectivity of the component 
(anything between i? sa 1 and i? ~ 0, while the anti-reflective coating is used to reduce the reflection 
at the second optical surface as much as possible so that this surface does not influence the light. 
Please note that the drawing is not to scale, the coatings are typically only a few microns thick on 
a several millimetre to centimetre thick substrate. 



the phase change upon reflection is either or 180°, depending on whether the second medium 
is optically thinner or thicker than the first. It is not shown here but Fresnel's equations can 
also be used to show that the phase change for the transmitted light at a lossless surface is zero. 
This contrasts with the definitions given in Section [2j] (see Figure ([3])ff.), where the phase shift 
upon any reflection is defined as zero and the transmitted light experiences a phase shift of tt/2. 
The following section explains the motivation for the latter definition having been adopted as the 
common notation for the analysis of modern optical systems. 

2.4.1 Composite optical surfaces 

Modern mirrors and beam splitters that make use of dielectric coatings are complex optical sys- 
tems, see Figure |8] whose reflectivity and transmission depend on the multiple interference inside 
the coating layers and thus on microscopic parameters. The phase change upon transmission or 
reflection depends on the details of the applied coating and is typically not known. In any case, the 
knowledge of an absolute value of a phase change is typically not of interest in laser interferometers 
because the absolute positions of the optical components are not known to sub-wavelength preci- 
sion. Instead the relative phase between the incoming and outgoing beams is of importance. In the 
following we demonstrate how constraints on these relative phases, i.e., the phase relation between 
the beams, can be derived from the fundamental principle of power conservation. To do this we 
consider a Michelson interferometer, as shown in Figure [9] with perfectly-reflecting mirrors. The 
beam splitter of the Michelson interferometer is the object under test. We assume that the magni- 
tude of the reflection r and transmission t are known. The phase changes upon transmission and 
reflection are unknown. Due to symmetry we can say that the phase change upon transmission ipt 
should be the same in both directions. However, the phase change on reflection might be different 
for either direction, thus, we write ipri for the reflection at the front and ipr2 for the reflection at 
the back of the beam splitter. 

Then the electric fields can be computed as 

Ei^r Eo e''^'-! ; E2 = t Eq e'^'. (18) 
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mirror 




mirror 



beam splitter 



E, 



Figure 9: The relation between the phase of the light field amplitudes at a beam splitter can be 
computed assuming a Michelson interferometer, with arbitrary arm length but perfectly-reflecting 
mirrors. The incoming field Eq is split into two fields Ei and E2 which are reflected atthe end 
mirrors and return to the beam splitter, as E3 and E4, to be recombined into two outgoing fields. 
These outgoing fields E5 and Eq are depicted by two arrows to highlight that these are the sum 
of the transmitted and reflected components of the returning fields. We can derive constraints 
for the phase of Ei and E2 with respect to the input field Eq from the conservation of energy: 



Eo' = E, 



Ee 



12 



We do not know the length of the interferometer arms. Thus, we introduce two further unknown 
phases: <i>i for the total phase accumulated by the field in the vertical arm and $2 for the total 
phase accumulated in the horizontal arm. The fields impinging on the beam splitter compute as 

Ey. = r Eo e'('^'-i+*i) ; E4 = t Eq e'('^'+*^^ (19) 

The outgoing fields are computed as the sums of the reflected and transmitted components: 

^ ' (20) 

Ee = £:o rt (e'(v't+v.i+*i) + gi (vt+v.2+*2)) ^ 

with i? = and T = t'^. 

It will be convenient to separate the phase factors into common and differential ones. We can 
write 

E5 = Eo e'"+ (R e'"- + T e"'"-) , (21) 



with 

and similarly 
with 



a+ = + + i ($1 + $2) ; a- = 95,-1 ~ipt + ^{'i>i- ^2) , (22) 
Eg = Eo rt e'^+ 2cos(^_), (23) 

l3+=ipt + ^ {(Prl + <^r2 + 'Si + ^'2) ; = \ ('^'^1 ~ "^''2 + *1 " ^2) ■ (24) 

For simplicity we now limit the discussion to a 50:50 beam splitter with r — t — l/-\/2, for which 
we can simplify the field expressions even further: 

E^ = Eo e'"+ cos(a_) ; E^ = Eo e''^+ cos(/3_). (25) 

Conservation of energy requires that |i?oP = l^'Sp + l^-eP, which in turn requires 

cos^(a_) + cos2 (/?_) = 1, (26) 

which is only true if 

a_-/3_ = (2iV + l)|, (27) 

with N as in integer (positive, negative or zero). This gives the following constraint on the phase 
factors 

-((prt+(p,2)-<^t = (2iV + l)|. (28) 

One can show that exactly the same condition results in the case of arbitrary (lossless) reflectivity 
of the beam splitter [IF . 

We can test whether two known examples fulflU this condition. If the beam-splitting surface 
is the front of a glass plate we know that ^pt — 0, ^pri = tt, (^^2 — 0, which conforms with 
Equation (28). A second example is the two-mirror resonator, see Section 2.2 If we consider the 
cavity as an optical 'black box', it also splits any incoming beam into a reflected and transmitted 
component, like a mirror or beam splitter. Further we know that a symmetric resonator must give 
the same results for fields injected from the left or from the right. Thus, the phase factors upon 
reflection must be equal Lpr = fri = fr2- The reflection and transmission coefficients are given by 
Equations ([t]) and Q as 

^ _ r2tjexp(-i2fci^) \ 

1 - rir2exp(-i2A:i:>)y ' ^ ' 
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and 

_ -tit2CXp{-ikD) 

1 — rir2 exp(— 1 2kD) 

We demonstrate a simple case by putting the cavity on resonance {kD — Nn). This yields 

r-cav = T"! - ; ^cav = Z , (31) 

with Tcav being purely real and ^cav imaginary and thus ipt = 7r/2 and ipr ^ which also agrees 



with Equation ( 28 ) 



In most cases we neither know nor care about the exact phase factors. Instead we can pick any 
set which fulfills Equation (|28|. For this document we have chosen to use phase factors equal to 
those of the cavity, i.e., ipt = 7r/2 and ipr — 0, which is why we write the reflection and transmission 
at a mirror or beam splitter as 

-Brofi Eo and E'trans = i t Eq. (32) 

In this definition r and t are positive real numbers satisfying — 1 for the lossless case. 

Please note that we only have the freedom to chose convenient phase factors when we do not 
know or do not care about the details of the optical system, which performs the beam splitting. 
If instead the details are important, for example when computing the properties of a thin coating 
layer, such as anti-reflex coatings, the proper phase factors for the respective interfaces must be 
computed and used. 



2.5 Lengths and tunings: numerical accuracy of distances 




■oi— 

A=lHm 



Figure 10: Illustration of an arm cavity of the Virgo gravitational- wave detector [56] : the macro- 
scopic length L of the cavity is approximately 3 km, while the wavelength of the NdiYAG laser 
is A ~ 1 /im. The resonance condition is only affected by the microscopic position of the wave 
nodes with respect to the mirror surfaces and not by the macroscopic length, i.e., displacement of 
one mirror by Ax = A/2 re-creates exactly the same condition. However, other parameters of the 
cavity, such as the flnesse, only depend on the macroscopic length L and not on the microscopic 
tuning. 



The resonance condition inside an optical cavity and the operating point of an interferometer 
depends on the optical path lengths modulo the laser wavelength, i.e., for light from an Nd:YAG 
laser length differences of less than 1 jim are of interest, not the full magnitude of the distances 
between optics. On the other hand, several parameters describing the general properties of an 
optical system, like the finesse or free spectral range of a cavity (see Section 5.11 depend on the 
macroscopic distance and do not change significantly when the distance is changed on the order of 
a wavelength. This illustrates that the distance between optical components might not be the best 
parameter to use for the analysis of optical systems. Furthermore, it turns out that in numerical 
algorithms the distance may suffer from rounding errors. Let us use the Virgo |56| arm cavities 
as an example to illustrate this. The cavity length is approximately 3 km, the wavelength is on 
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the order of 1 jini, the mirror positions are actively controlled with a precision of 1 pm and the 
detector sensitivity can be as good as 10"^* m, measured on ^ 10 nis timescales (i.e., many samples 
of the data acquisition rate) . The floating point accuracy of common, fast numerical algorithms 
is typically not better than 10"^^. If we were to store the distance between the cavity mirrors as 
such a floating point number, the accuracy would be limited to 3 pm, which does not even cover 
the accuracy of the control systems, let alone the sensitivity. 

A simple and elegant solution to this problem is to split a distance D between two optical 
components into two parameters [30]: one is the macroscopic 'length' L, defined as the multiple 
of a constant wavelength Aq yielding the smallest difference to D. The second parameter is the 
microscopic tuning T that is defined as the remaining difference between L and D, i.e., D = L + T. 
Typically, Xq can be understood as the wavelength of the laser in vacuum, however, if the laser 
frequency changes during the experiment or multiple light fields with different frequencies are used 
simultaneously, a default constant wavelength must be chosen arbitrarily. Please note that usually 
the term A in any equation refers to the actual wavelength at the respective location as X = Xo/n 
with n the index of refraction at the local medium. 

We have seen in Section [2. 1| that distances appear in the expressions for electromagnetic waves 
in connection with the wave number, for example, 

E2 = El exp{~ikz). (33) 

Thus, the difference in phase between the field at z — zi and z ~ zi + D is given as 

(/3 = -kD. (34) 

We recall that k = 2tt/X = uj/c. We can define ljq — 27r c/Aq and k^ = ojq/c. For any given 
wavelength A we can write the corresponding frequency as a sum of the default frequency and a 



difference frequency w = Wq + Aw. Using these definitions, we can rewrite Equation (34) with 
length and tuning as 

, ^ luqL AluL luqT AluT 

-ip = kD = — + + — + . 35 

c c c c 

The first term of the sum is always a multiple of 2tt, which is equivalent to zero. The last term 
of the sum is the smallest, approximately of the order Au ■ 10^^^. For typical values of L « 1 m, 
T <1 iim and Aw < 27r • 100 MHz we find that 

^^0, ^<2, ^<6, ^<210-^ (36) 
c c ~ c ~ c 

which shows that the last term can often be ignored. 

We can also write the tuning directly as a phase. We define as the dimensionless tuning 

= iOoT/c. (37) 

This yields 



(i-T)=exp i^T^Uexp i^0 . (38) 



exp 

' " ' ' ' c uqJ " V Wo 
The tuning <p is given in radian with 27r referring to a microscopic distance of one wavelengttj^ Aq . 

^Note that in other pubUcations the tuning or equivalent microscopic displacements are sometimes defined via 
an optical path-length difference. In that case, a tuning of 2it is used to refer to the change of the optical path 
length of one wavelength, which, for example, if the reflection at a mirror is described, corresponds to a change of 
the mirror's position of Ao/2. 
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Finally, we can write the following expression for the phase difference between the light field 
taken at the end points of a distance D: 



^ = .kD = -[^ + ^^], (39) 

c Wo 



or if we neglect the last term from Equation (36) we can approximate (w/wq ~ 1) to obtain 



(40) 

This convention provides two parameters L and (f), that can describe distances with a markedly 
improved numerical accuracy. In addition, this definition often allows simplification of the algebraic 
notation of interferometer signals. By convention we associate a length L with the propagation 
through free space, whereas the tuning will be treated as a parameter of the optical components. 
Effectively the tuning then represents a microscopic displacement of the respective component. If, 
for example, a cavity is to be resonant to the laser light, the tunings of the mirrors have to be the 
same whereas the length of the space in between can be arbitrary. 
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2.6 Revised coupling matrices for space and mirrors 

Using the definitions for length and tunings we can rewrite the couphng equations for mirrors 



and spaces introduced in Section 2.1 as follows. The mirror coupling becomes (compare this to 



"4 



\t rexp(i2(jl)- 
rexp(-i2^i^) it 



"3 
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reference plane 



as 



Figure 11: 

Figure |6]), and the amplitude coupling for a 'space', formally written as in Figure [t] is now written 



as 



space (L, n) 



a2 

04 



exp(— iAfcnL) \ f '^^ 

exp(— i AfcnL) J \ 



ai 
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Figure 12: 
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2.7 Finesse examples 

2.7.1 Mirror reflectivity and transmittance 



fexample_rt_mirror 
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Figure 13: Finesse example: Mirror reflectivity and transmittance. 

We use Finesse to plot the amplitudes of the light fields transmitted and reflected by a mirror 
(given by a single surface). Initially, the mirror has a power reflectance and transmittance of 
R — T — 0.5 and is, thus, lossless. For the plot in Figure [T3| we tune the transmittance from 0.5 to 
0. Since we do not explicitly change the reflectivity, R remains at 0.5 and the mirror loss increases 
instead, which is shown by the trace labelled 'total' corresponding to the sum of the reflected and 
transmitted light power. The plot also shows the phase convention of a 90° phase shift for the 
transmitted light. 

Finesse input file for 'Mirror reflectivity and transmittance' 

laser 11 1 nl 7, laser with P=1W at the default frequency 
space si 1 nl n2 space of Im length 

mirror ml 0.5 0.5 n2 nS 7, mirror with T=R=0 . 5 at zero tuning 
ad ad_t nS 7« an 'amplitude' detector for transmitted light 

ad ad_r n2 7o an 'amplitude' detector for reflected light 

set t ad_t abs 
set r ad_r abs 

func total = $r~2 + $t~2 7« computing the sum of the reflected and transmitted power 



xaxis ml t lin 0.5 100 7« changing the transmittance of the mirror 'ml' 
yaxis abstdeg 7. plotting amplitude and phase of the results 
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Figure 14: Finesse example: Length and tunings. 



2.7.2 Length and tunings 

This Finesse file demonstrates the conventions for lengths and microscopic positions introduced 
in Section |2.5| The top trace in Figure [M] depicts the phase change of a beam reflected by a 
beam splitter as the function of the beam splitter tuning. By changing the tuning from to 180° 
the beam splitter is moved forward and shortens the path length by one wavelength, which by 
convention increases the light phase by 360°. On the other hand, if a length of a space is changed, 
the phase of the transmitted light is unchanged (for the default wavelength Afc = 0), as shown the 
in the lower trace. 

Finesse input file for 'Length and tunings' 

laser 11 1 nl 7, laser with P=1W at the default frequency 
space si 1 1 nl n2 */, space of Im length 

bs bl 1 n2 n3 dump dump '/, beam splitter as 'turning mirror' , normal incidence 

space s2 1 1 nS n4 7, another space of Im length 
ad adl n4 7« amplitude detector 

7. for the plot we perform two sequenctial runs of Finesse using 'mkat' 

7. 1) first trace: change microscopic position of beamsplitter 

runl: xaxis bl phi lin 180 100 

7 2) second trace: change length of space si 

run2: xaxis si L lin 1 2 100 

yaxis deg 7. plotting the phase of the results 
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3 Light with Multiple Frequency Components 



So far we have considered the electromagnetic field to be monochromatic. This has allowed us to 
compute light-field amplitudes in a quasi-static optical setup. In this section, we introduce the 
frequency of the light as a new degree of freedom. In fact, we consider a field consisting of a finite 
and discrete number of frequency components. We write this as 



Eit,z) 



E 

3 



aj exp (i {ujj t — kjz)), 



(41) 



with complex amplitude factors aj, uij as the angular frequency of the light field and kj = ujj/c. 
In many cases the analysis compares different fields at one specific location only, in which case we 
can set z = and write 

£;(t) = flj exp(iwji). (42) 



In the following sections the concept of light modulation is introduced. As this inherently involves 
light fields with multiple frequency components, it makes use of this type of field description. Again 
we start with the two-mirror cavity to illustrate how the concept of modulation can be used to 
model the effect of mirror motion. 



a) 



m 



b) 



/(*) 



f(t) = sin(a;t -|- msin(fli)) 

> t 



f[t) = [l-l-msin(nt)] sin(tjt) 

r 

> t 



Figure 15: Example traces for phase and amplitude modulation: the upper plot a) shows a phase- 
modulated sine wave and the lower plot b) depicts an amplitude-modulated sine wave. Phase 
modulation is characterised by the fact that it mostly affects the zero crossings of the sine wave. 
Amplitude modulation affects mostly the maximum amplitude of the wave. The equations show 
the modulation terms in red with m the modulation index and the modulation frequency. 



3.1 Modulation of light fields 

Laser interferometers typically use three different types of light fields: the laser with a frequency of, 
for example, / « 2.8- 10^"* Hz, radio frequency (RF) sidebands used for interferometer control with 
frequencies (offset to the laser frequency) of / « 1 • 10^ to 150 • 10® Hz, and the signal sidebands 
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at frequencies of 1 to 10,000 As these modulations usually have as their origin a change in 

optical path length, they are often phase modulations of the laser frequency, the RF sidebands are 
utilised for optical readout purposes, while the signal sidebands carry the signal to be measured 
(the gravitational- wave signal plus noise created in the interferometer). 

Figure [Ts] shows a time domain representation of an electromagnetic wave of frequency ujq, 
whose amplitude or phase is modulated at a frequency Q. One can easily see some characteristics 
of these two types of modulation, for example, that amplitude modulation leaves the zero crossing 
of the wave unchanged whereas with phase modulation the maximum and minimum amplitude of 
the wave remains the same. In the frequency domain in which a modulated field is expanded into 
several unmodulated field components, the interpretation of modulation becomes even easier: any 
sinusoidal modulation of amplitude or phase generates new field components, which are shifted 
in frequency with respect to the initial field. Basically, light power is shifted from one frequency 
component, the carrier, to several others, the sidebands. The relative amplitudes and phases of 
these sidebands differ for different types of modulation and different modulation strengths. This 
section demonstrates how to compute the sideband components for amplitude, phase and frequency 
modulation. 



3.2 Phase modulation 

Phase modulation can create a large number of sidebands. The number of sidebands with noticeable 
power depends on the modulation strength (or depth) given by the modulation index m. Assuming 
an input field 

Ein = Eo exp {i UJQ t), (43) 
a sinusoidal phase modulation of the field can be described as 

E ^ Eq exp(^i{ujot + mcos{^lt))y (44) 

This equation can be expanded using the identity [25] 

oo 

exp{i z COS ip) — ^ i'' Jk{z) exp{i kip), (45) 

fc — — oo 

with Bessel functions of the first kind Jk{rn). We can write 

oo 

E = Eq exp(iwoi) ^ Jk{m) exp{ikU.t). (46) 

k— — oo 

The field for fc = 0, oscillating with the frequency of the input field ujq, represents the carrier. The 
sidebands can be divided into upper {k > 0) and lower (k < 0) sidebands. These sidebands are 
light fields that have been shifted in frequency by k fi. The upper and lower sidebands with the 
same absolute value of k are called a pair of sidebands of order k. Equation (46 1 shows that the 
carrier is surrounded by an infinite number of sidebands. However, for small modulation indices 
(m < 1) the Bessel functions rapidly decrease with increasing k (the lowest orders of the Bessel 
functions are shown in Figure 16). For small modulation indices we can use the approximation |5] 

*('")^(f)'E;fe^^^(?)' + o('"'«)- (-) 



^The signal sidebands are sometimes also called audio sidebands because of their frequency range. 
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In which case, only a few sidebands have to be taken into account. For m ^ 1 we can write 



E — Eq exp (i luq t) 

X ^Jo(to) — i J_i(m) exp (— i t) + i Ji(to) exp(if2f)^, 



(48) 



and with 

J_fe(m) = Jk{m), (49) 

we obtain 

E ^ Eq exp{iLUot) (^1 + i ^ (^exp (-i t) + exp (i t)^ ^ , (50) 

as the first-order approximation in m. In the above equation the carrier field remains unchanged 
by the modulation, therefore this approximation is not the most intuitive. It is clearer if the 
approximation up to the second order in m is given: 

/ 771"^ 771 / \ \ 

E = Eq eTip{iujot) fl — +i — (^e^p{-int) +e^p{int)j j , (51) 

which shows that power is transferred from the carrier to the sideband fields. 

Higher-order expansions in m can be performed simply by specifying the highest order of Bessel 
function, which is to be used in the sum in Equation (46), i.e., 

order 

E ^ Eo exp (iojat) ^ i'' Jk (m) exp (ikflt). (52) 

k— — order 




Figure 16: Some of the lowest-order Bessel functions Jk{x) of the first kind. For small x the 
expansion shows a simple x'^ dependency and higher-order functions can often be neglected. 



3.3 Frequency modulation 

For small modulation, indices, phase modulation and frequency modulation can be understood as 
different descriptions of the same effect [30^ . Following the same spirit as above we would assume 
a modulated frequency to be given by 

oj — ujQ + m' cos (fit), (53) 
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and then we might be tempted to write 



E ^ Eq exp (^1 {ujQ + m' cos {Q t))t^, (54) 
which would be wrong. The frequency of a wave is actually defined as lu/{2tt) ~ f — dtp/dt. Thus, 



to obtain the frequency given in Equation (53), we need to have a phase of 



771 

uJot+ sin (fit). (55) 
For consistency with the notation for phase modulation, we define the modulation index to be 

with Aw as the frequency swing - how Jar the frequency is shifted by the modulation - and SI the 
modulation frequency - how fast the frequency is shifted. Thus, a sinusoidal frequency modulation 
can be written as 

E = Eo exp (i If) — Eq exp M fujo t + cos {Q t)j j , (57) 



which is exactly the same expression as Equation (44 1 for phase modulation. The practical dif- 
ference is the typical size of the modulation index, with phase modulation having a modulation 
index of m < 10, while for frequency modulation, typical numbers might be to > 10'*. Thus, in the 
case of frequency modulation, the approximations for small to are not valid. The series expansion 



using Bessel functions, as in Equation (46), can still be performed, however, very many terms of 



the resulting sum need to be taken into account. 
3.4 Amplitude modulation 

In contrast to phase modulation, (sinusoidal) amplitude modulation always generates exactly two 
sidebands. Furthermore, a natural maximum modulation index exists: the modulation index is 
defined to be one (to = 1) when the amplitude is modulated between zero and the amplitude of 
the unmodulated field. 

If the amplitude modulation is performed by an active element, for example by modulating the 
current of a laser diode, the following equation can be used to describe the output field: 



E ^ E, 



exp (i wq i) ^1 + TO cos (f2 1) 
Eq exp (i Wo t) exp (— i n t) 



(58) 



However, passive amplitude modulators (like acousto-optic modulators or electro-optic modulators 
with polarisers) can only reduce the amplitude. In these cases, the following equation is more 
useful: 

E ^ £^0 exp (i Wot) f 1 - f f 1 - cos (fi t)^ 



(59) 

= Eo cxp(iwo<) (1-f + f cxp{int) + f exp(-if7i) 
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Figure 18: Amplitude and phase modulation in the 'phasor' picture. The upper plots a) illustrate 
how a phasor diagram can be used to describe phase modulation, while the lower plots b) do the 
same for amplitude modulation. In both cases the left hand plot shows the carrier in blue and 
the modulation sidebands in green as snapshots at certain time intervals. One can see clearly 
that the upper sideband (cjq + ^) rotates faster than the carrier, while the lower sideband rotates 
slower. The right plot in both cases shows how the total field vector at any given time can be 
constructed by adding the three field vectors of the carrier and sidebands. [Drawing courtesy of 
Simon Chelkowski] 
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3.5 Sidebands as phasors in a rotating frame 



A common method of visualising the behaviour of sideband fields in interferometers is to use phase 
diagrams in which each field amplitude is represented by an arrow in the complex plane. 

We can think of the electric field amplitude Eoexp(iujQt) as a vector in the complex plane, 
rotating around the origin with angular velocity ljq. To illustrate or to help visualise the addition 
of several light fields it can be useful to look at this problem using a rotating reference frame, 
defined as follows. A complex number shall be defined as z = x + iy so that the real part is plotted 
along the x-axis, while the y-axis is used for the imaginary part. We want to construct a new 
coordinate system {x' , y') in which the field vector is at a constant position. This can be achieved 
by defining 

X ~ x' cosioot^y' smujQt 
y — x' smujQt + y' cosujQt, 



(60) 



or 



x' =x cos(— Wot) — 2/ sin(— wq^) 
y' = X sin (— cjqO + y cos (— Wq^) . 



(61) 



Figure 17 illustrates how the transition into the rotating frame makes the field vector to appear 
stationary. The angle of the field vector in a rotating frame depicts the phase offset of the field. 
Therefore these vectors are also called phasors and the illustrations using phasors are called phasor 
diagrams. Two more complex examples of how phasor diagrams can be employed is shown in 
Figure [18] fTT. 

Phasor diagrams can be especially useful to see how frequency coupling of light field ampli- 
tudes can change the type of modulation, for example, to turn phase modulation into amplitude 
modulation. An extensive introduction to this type of phasor diagram can be found in |40j . 



mirror 
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Figure 19: A sinusoidal signal with amplitude Os frequency ujs and phase offset ips is applied to 
a mirror position, or to be precise, to the mirror tuning. The equation given for the tuning (j) 



assumes that uJs/ojo ^ 1, see Section 2.5 



3.6 Phase modulation through a moving mirror 

Several optical components can modulate transmitted or reflected light fields. In this section we 
discuss in detail the example of phase modulation by a moving mirror. Mirror motion does not 
change the transmitted light; however, the phase of the reflected light will be changed as shown in 
Equation (11 1. 

We assume sinusoidal change of the mirror's tuning as shown in Figure |19[ The position 
modulation is given as x^^ = Og cos{uJst + (ps), and thus the reflected field at the mirror becomes 
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(assuming 04 ~ 0) 

03 = r fli exp(— i 200) exp (i 2fca;ni) ~ Ta^ exp(— i 2(/)o) exp 2kQas^ cos^ujgt + ips)j , (62) 
setting m = 2/0903. This can be expressed as 

03 = raiexp(-i2(/)o) (l + if exp(^-i(a;st + <y5s)^ +if^exp(^i(wst + (/3s)^^ 

= raiexp(-i20o) (1 + f exp (^-i(w3i + v?s - 7r/2)) (63) 
+ f exp (^i {ujst + ips + 7r/2)j j . 

3.7 Coupling matrices for beams with multiple frequency components 

The coupUng between electromagnetic fields at optical components introduced in Section [2] referred 
only to the amplitude and phase of a simplified monochromatic field, ignoring all the other param- 
eters of the electric field of the beam given in Equation ([!]). However, this mathematical concept 
can be extended to include other parameters provided that we can find a way to describe the total 
electric field as a sum of components, each of which is characterised by a discrete value of the 
related parameters. In the case of the frequency of the light field, this means we have to describe 
the field as a sum of monochromatic components. In the previous sections we have shown how this 
could be done in the special case of an initial monochromatic field that is subject to modulation: 
if the modulation index is small enough we can limit the amount of frequency components that 
we need to consider. In many cases it is actually sufficient to describe a modulation only by the 
interaction of the carrier at ujq (the unmodulated field) and two sidebands with a frequency offset 
of ±Wm to the carrier. A beam given by the sum of three such components can be described by a 
complex vector: 

a('^o) \ / flf^o \ 
a{uJo - LU,n) = a^i (64) 
a{uJo +UJ,n) / \ a^2 / 

Lul and ujQ + ujm — w2. In the case of a phase modulator that applies a 
modulation of small modulation index m to an incoming light field di , we can describe the coupling 
of the frequency component as follows: 



with ujQ = wO, ujQ 



which can be written in matrix form: 



Jo{m) Ji{m) J_i(m) 
d2 = I J-i{m) Ja{m) 

Ji (m) Jo [m) 



ai. 



(65) 



(66) 



And similarly, we can write the complete coupling matrix for the modulator component, for ex- 
ample, as 



/ 0,2,wVl \ 

a2,u)i 

14,11)0 

a4,toi 

V a4,tu2 / 



J_i(m) 
Ji(m) 



V 



Ji(m) 
Joim) 











J_i(m) 


Joim) 










Jo(m) 
J_i(m) 
Ji(m) 







Ji(to) 
Jn{m) 




\ 




J_i(to) 


Jo(m) J 



( Oi^^o \ 
0-\,w2 

Mil 

\ 03,1112 / 



(67) 
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3.8 Finesse examples 
3.8.1 Modulation index 



fexample_bessel 
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Figure 20: Finesse example: Modulation index. 



This file demonstrates the use of a modulator. Phase modulation (with up to five higher harmonics 
is applied to a laser beam and amplitude detectors are used to measure the field at the first three 



harmonics. Compare this to Figure 16 as well. 
Finesse input file for 'Modulation index' 



laser il 1 nO 7, laser P=1W f_offset=OHz 

mod eoml 40k .05 5 pm nO nl phase modulator f_mod=40kHz, modulation index=0.05 

ad bessell 40k nl 7, amplitude detector f=40kHz 

ad bessel2 80k nl 7, amplitude detector f=80kHz 

ad besselS 120k nl '/, amplitude detector f=120kHz 

xaxis eoml midx lin 10 1000 x-axis: modulation index of eoml 



yaxis abs 



7 y-axis: plot 'absolute' amplitude 



3.8.2 Mirror modulation 

Finesse offers two different types of modulators: the 'modulator' component shown in the exam- 
ple above, and the 'fsig' command, which can be used to apply a signal modulation to existing 
optical components. The main difference is that 'fsig' is meant to be used for transfer function 
computations. Consequently Finesse discards all nonlinear terms, which means that the sideband 
amplitude is proportional to the signal amplitude and harmonics are not created. 

Finesse input file for 'Mirror modulation' 

laser il 1 nl 1 laser P=1W f_offset=OHz 

space si 1 1 nl n2 °/t space of Im length 
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fexample_mirror_modulation Mon Aug 17 10:35:43 2009 
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Figure 21: Finesse example: Mirror modulation. 

bs bl 1 n2 n3 dump dump 7. beam splitter as 'turning mirror' , normal incidence 

space s2 1 1 n3 n4 '/, another space of Im length 

fsig sigl bl 40k 1 7, signal modulation applied to beam splitter bl 

ad upper 40k n4 7, amplitude detector f=40kHz 

ad lower -40k n4 7, amplitude detector f=-40kHz 

ad harmonic 80k n4 7, amplitude detector f=80kHz 

xaxis sigl amp lin 1 10 100 7o x-axis: amplitude of signal modulation 

yaxis abs 7. y-axis: plot 'absolute' amplitude 
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4 Optical Readout 



In previous sections we have dealt with the amphtude of hght fields directly and also used the 
amplitude detector in the FINESSE examples. This is the advantage of a mathematical analysis 
versus experimental tests, in which only light intensity or light power can be measured directly. 
This section gives the mathematical details for modelling photo detectors. 

The intensity of a field impinging on a photo detector is given as the magnitude of the Poynting 
vector, with the Poynting vector given as [35] 

S = Ex H = —E X B. (68) 
Mo 

Inserting the electric and magnetic components of a plane wave, we obtain 

= —E^ = ceoEl cos^{ujt) = '^E^ (1 + cos(2a;i)) , (69) 
Moc 2 

with £0 the electric permeability of vacuum and c the speed of light. 

The response of a photo detector is given by the total flux of effective radiatiorj^ during the 
response time of the detector. For example, in a photo diode a photon will release a charge in the 
n-p junction. The response time is given by the time it takes for the charge to travel through the 
detector (and further time may be taken up in the electronic processing of the signal). The size of 
the photodiode and the applied bias voltage determine the travel time of the charges with typical 
values of approximately 10 ns. Thus, frequency components faster than perhaps 100 MHz are not 
resolved by a standard photodiode. For example, a laser beam with a wavelength of A 1064 nm 
has a frequency of / = c/A « 282 10^^ Hz = 282 THz. Thus, the 2a; component is much too fast 
for the photo detector; instead, it returns the average power 

1^1 = '^El (70) 

In complex notation we can write 

\§\ ^ ^EE*. (71) 

However, for more intuitive results the light fields can be given in converted units, so that the 
light power can be computed as the square of the light field amplitudes. Unless otherwise noted, 
throughout this work the unit of light field amplitudes is y^watt. Thus, the notation used in this 
document to describe the computation of the light power of a laser beam is 

P^EE*. (72) 



4.1 Detection of optical beats 

What is usually called an optical beat or simply a beat is the sinusoidal behaviour of the intensity 
of two overlapping and coherent fields. For example, if we superpose two fields of slightly different 
frequency, we obtain 

E = Eq cos{uJit) + Eq cos{uj2t) 

P = E"^ ^El (cos2(a;ii) + cos2(w20 + 2cos(a;it) cos(w2t)) ^'^^^ 
— Eq (cos^(a;ii) + cos^(w2t) + cos(w+t) + cos(w_t)) , 

*Thc term effective refers to that amount of incident light, which is converted into photo-electrons that are then 
usefully extracted from the junction (i.e., do not recombine within the device). This fraction is usually referred to 
as quantum efficiency rj of the photo diode. 
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with w+ = wi + W2 and w_ — uji — uj2- In this equation the frequency w„ can be very small and 



can then be detected with the photodiode as illustrated in Figure 22 



Pdiodo -i^o (l + cos(w_t)) (74) 

Using the same example photodiode as before: in order to be able to detect an optical beat a;_ 
would need to be smaller than 100 MHz. If we take two, sightly detuned Nd:YAG lasers with 
/ — 282 THz, this means that the relative detuning of these lasers must be smaller than 10"''. 
In general, for a field with several frequency components, the photodiode signal can be written 

as 

N N 

\E\^ = = ^ ^ a,a* e'("'-"^)*. (75) 

i=0 j=0 

For example, if the photodiode signal is filtered with a low-pass filter, such that only the DC part 
remains, we can compute the resulting signal by looking for all components without frequency 
dependence. The frequency dependence vanishes when the frequency becomes zero, i.e., in all 
parts of Equation ( |75[ ) with uji ~ ojj. The output is a real number, calculated like this: 

x^^^^^a^a* with {z, j | i,j e {0, . . . , A^} A u!i=ujj}. (76) 



4.2 Signal demodulation 

A typical application of light modulation, is its use in a modulation-demodulation scheme, which 
applies an electronic demodulation to a photodiode signal. A 'demodulation' of a photodiode signal 
at a user-defined frequency lo^, performed by an electronic mixer and a low-pass filter, produces a 
signal, which is proportional to the amplitude of the photo current at DC and at the frequency 
ujo ± uJx- Interestingly, by using two mixers with different phase offsets one can also reconstruct 
the phase of the signal, or to be precise the phase difference of the light at ujq ± with respect to 
the carrier light. This feature can be very powerful for generating interferometer control signals. 

Mathematically, the demodulation process can be described by a multiplication of the output 
with a cosine: cos{u;x + ipx) {^x is the demodulation phase), which is also called the 'local oscillator'. 
After the multiplication was performed only the DC part of the result is taken into account. The 
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signal is 

N N 

So = I^P = =^^a,a* e't"''"^'*. (77) 

i=0 j=0 



(78) 



Multiplied with the local oscillator it becomes 

51 = S-o •cos(a;:,i + ¥':.) = (e'("-*+'^-) +e-'("-*+'^-)) 

N N 
z=0 j=0 

With = aia* and e''^'^ ' = gi(w,-a;j)t ^^g^^^ ^^.j^.^ 

/ N N N \ 

^2=0 i=0 j=i+l J 



When looking for the DC components of Si we get the following 121]: 

'5'i,DC = Y.\{A,, e-'-^- + A*^ e'v^) with {z, j j e {0, . . . , iV} A c^,, = ojJ 



(80) 



This would be the output of a mixer and a subsequent low-pass filter. The results for ifx~() and 
= 7r/2 are called in-phase and in- quadrature, respectively (or also first and second quadrature). 
They are given by 

'S'l, DC, phase — Yl'^i'^ij} 1 

S"!, DC, quad — E^l^ij}- ^ ^ 

If only one mixer is used, the output is always real and is determined by the demodulation phase. 
However, with two mixers generating the in-phase and in-quadrature signals, it is possible to 
construct a complex number representing the signal amplitude and phase: 

z = ^a,a* with {iJ\i,j€{0,...,N}Auj,j=uJx}. (82) 

Often several sequential demodulations are applied in order to measure very specific phase in- 
formation. For example, a double demodulation can be described as two sequential multiplications 
of the signal with two local oscillators and taking the DC component of the result. First looking 
at the whole signal, we can write: 

'5'2 = 5*0 • cos(a;a;i + ip^) cos{ujyt + ipy). (83) 

This can be written as 

<S'2 = So^{cOs{ujyt + LOxt + ify + ipx) + COs{ujyt - Uxt + ify - ipx)) .„.^ 

= S'o^(cos(a;+t + (^+) + cos(a;_t + (/J_)), 

and thus reduced to two single demodulations. Since we now only care for the DC component we 
can use the expression from above (Equation (82 1). These two demodulations give two complex 
numbers: 

zl = E^ij with {i,j\i,je{0,...,N} A uji ~ ujj = uj+}, 

z2 = J^Au with {k,l\ k,le{0,...,N} A ujk-uji=uj^}. ^^^^ 
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The demodulation phases are apphed as fohows to get a real output (two sequential mixers) 

x = ^{{zie-'"^''+Z2e''^-)e-'"^y}. (86) 

In a typical setup, a \iscr-dcfincd demodulation phase for the first frequency (here ipx) is given. If 
two mixers are used for the second demodulation, we can reconstruct the complex number 

z = zi e-''^'' +Z2 e''^''. (87) 

More demodulations can also be reduced to single demodulations as above. 
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4.3 Finesse examples 
4.3.1 Optical beat 



fexample_optical_beat 



Mon Aug 17 21:02:00 2009 



■ pd1 n6 ' 

■ ad10k n6 
-ad1 n6 , 
-pd10k n6 
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Figure 23: Finesse example: Optical beat. 

In this example two laser beams are superimposed at a 50:50 beam splitter. The beams have a 
slightly different frequency: the second beam has a 10 kHz offset with respect to the first (and to 
the default laser frequency). The plot illustrates the output of four different detectors in one of 
the beam splitter output ports, while the phase of the second beam is tuned from 0° to 180°. The 
photodiode 'pdl' shows the total power remaining constant at 1. The amplitude detectors 'adl' 
and 'adlOk' detect the laser light at Hz (default frequency) and 10 kHz respectively. Both show 
a constant absolute of a/1/2 and the detector 'adlOk' tracks the tuning of the phase of the second 
laser beam. Finally, the detector 'pdlOk' resembles a photodiode with demodulation at 10 kHz. In 
fact, this represents a photodiode and two mixers used to reconstruct a complex number as shown 
in Equation ( 82 1 . One can see that the phase of the resulting electronic signal also directly follows 



the phase difference between the two laser beams. 
Finesse input file for 'Optical beat' 



const freq 10k creating a constant for the frequency offset 

laser 11 1 nl laser with P=1W at the default frequency 

space si In 1 nl n2 '/, space of Inm length 

laser 12 1 $freq n3 7. a second laser with f=10kHz frequency offset 

space s2 In 1 n3 n4 7, another space of Inm length 

bs bl 0.5 0.5 n2 n5 dump n4 "/. 50:50 beam splitter 

space s3 In 1 n5 n6 another space of Inm length 

ad adO n6 '/, amplitude detector at f=0Hz 

ad adlOk $freq n6 '/« amplitude detector at f=10kHz 

pd pdl n6 y, simple photo detector 

pdl pdlOk $freq n6 Z photo detector with demodulation at lOkHz 
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xaxis 12 phi lin 180 100 changing the phase of the 12-beam 
yaxis abstdeg plotting amplitude and phase 
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5 Basic Interferometers 



The large interferometric gravitational-wave detectors currently in operation are based on two fun- 
damental interferometer topologies: the Fabry-Perot and the Michelson interferometer. The main 
instrument is very similar to the original interferometer concept used in the famous experiment 
by Michelson and Morley, published in 1887 |43]. The main difference is that modern instru- 
ments use laser light to illuminate the interferometer to achieve much higher accuracy. Already 
the first prototype by Forward and Weiss has thus achieved a sensitivity a million times better 
than Michelson's original instrument '19]. In addition, in current gravitational- wave detectors, the 
Michelson interferometer has been enhanced by resonant cavities, which in turn have been derived 
from the original idea for a spectroscopy standard published by Fabry and Perot in 1899 [T7]- The 
following section will describe the fundamental properties of the Fabry-Perot interferometer and 
the Michelson interferometer. A thorough understanding of these basic instruments is essential for 
the study of the high-precision interferometers used for gravitational-wave detection. 
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Figure 24: Typical optical layout of a two-mirror cavity, also called a Fabry-Perot interferometer. 
Two mirrors form the Fabry-Perot interferometer, a laser beam is injected through one of the 
mirrors and the reflected and transmitted light can be detected by photo detectors. 



5.1 The two-mirror cavity: a Fabry-Perot interferometer 

We have computed the field amplitudes in a linear two-mirror cavity, also called Fabry-Perot 
interferometer, in Section [2. 2[ In order to understand the features of this optical instrument it is 
of interest to have a closer look at the power circulation in the cavity. A typical optical layout is 



shown in Figure 24 two parallel mirrors form the Fabry-Perot cavity. A laser beam is injected 
through the first mirror (at normal incidence). 

The behaviour of the (ideal) cavity is determined by the length of the cavity L, the wavelength 
of the laser A and the reflectivity and transmittance of the mirrors. Assuming an input power of 
|aop = 1, we obtain 



Pi = |ai| 



1 + R1R2 - 2rir2 cos {2kL) '' 



with k = 27r/A, P, T = t and R = r , as defined in Section 1.4 Similarly we could compute the 

nple, 

(89) 



transmission of the optical system as the input-output ratio of the field amplitudes. For example, 

02 — ii<2 cxp(— i fcL) 



ao 1 — rir2 exp(— i 2fcL) 

is the frequency-dependent transfer function of the cavity in transmission (the frequency depen- 
dency is hidden inside the k = 27r//c). 



Figure 25 shows a plot of the circulating light power Pi over the laser frequency. The maximum 



power is reached when the cosine function in the denominator becomes equal to one, i.e., at 
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Figure 25: Power enhancement in a two-mirror cavity as a function of the laser-hght frequency. 
The peaks marks the resonances of the cavity, i.e., modes of operation in which the injected hght 
is resonantly enhanced. The frequency distance between two peaks is called free-spectral range 
(FSR). 



kL — Ntt with N an integer. This is called the cavity resonance. The lowest power values are 
reached at anti-resonance when kL — {N + l/2)7r. We can also rewrite 

2fcL = .^^2./^.||(, (90) 
c c i^ bK 



with FSR being the free-spectral range of the cavity as shown in Figure 25 Thus, it becomes clear 
that resonance is reached for laser frequencies 

fr^N- FSR, (91) 

where N is an integer. 

Another characteristic parameter of a cavity is its linewidth, usually given as full width at half 
maximum (FWHM) or its pole frequency, fp. In order to compute the linewidth we have to ask at 
which frequency the circulating power becomes half the maximum: 

This results in the following expression for the full linewidth: 

FWHM ^2fp^?^ arcsin ( . (93) 

The ratio of the linewidth and the free spectral range is called the finesse of a cavity: 

F^^^ = ^ (94) 

FWHM 2 arcsin fl-SmV ^ ^ 
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Figure 26: This figure compares the fields reflected by, transmitted by and circulating in a Fabry- 
Perot cavity for the three different cases; over-coupled, under-coupled and impedance matched 
cavity (in all cases Ti +T2 = 0.2 and the round-trip loss is 1%). The traces show the phase and 
amplitude of the electric field as a function of laser frequency detuning. 
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In the case of high finesse, i.e., ri and r2 are close to 1 we can use the fact that the argument of 
the arcsin function is small and make the approximation 



F^-^^^^- . (95) 

1 - rir2 1 - rir2 
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Figure 27: Power transmitted and circulating in a two mirror cavity with input power 1 W. The 
mirror transmissions are set such that Ti + T2 = 0.8 and the reflectivities of both mirrors are set 
as i? = 1 — T. The cavity is undercoupled for Ti < 0.4, impedance matched at Ti = T2 = 0.4 and 
overcoupled for Ti > 0.4. The transmission is maximised in the impedance-matched case and falls 
similarly for over or undercoupled settings. However, the circulating power (and any resonance 
performance of the cavity) is much larger in the overcoupled case. 

The behaviour of a two mirror cavity depends on the length of the cavity (with respect to the 
frequency of the laser) and on the reflectivities of the mirrors. Regarding the mirror parameters 
one distinguishes three casetj^ 

• when Ti < T2 the cavity is called undercoupled 

• when Ti — T2 the cavity is called impedance matched 

• when Ti > T2 the cavity is called overcoupled 

The differences between these three cases can seem subtle mathematically but have a strong 
impact on the application of cavities in laser systems. One of the main differences is the phase 
evolution of the light flelds, which is shown in Figure |26| The circulating power shows that the 
resonance effect is better used in over-coupled cavities; this is illustrated in Figure [27j which 
shows the transmitted and circulating power for the three different cases. Only in the impedance- 
matched case can the cavity transmit (on resonance) all the incident power. Given the same total 



^Please note that in the presence of losses the coupling is defined with respect to the transmission and losses. In 
particular, the impedance-matched case is defined as Ti = T2 ■ Loss, so that the input power transmission exactly 
matches the light power lost in one round-trip. 
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transmission Ti + T2, the overcoupled case allows for the largest circulating power and thus a 
stronger 'resonance effect' of the cavity, for example, when the cavity is used as a mode filter. 
Hence, most commonly used cavities are impedance matched or overcoupled. 
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Figure 28: Typical optical layout of a Michelson interferometer: a laser beam is split into two 
and sent along two perpendicular interferometer arms. We will label the directions in a Michelson 
interferometer as North, East, West and South in the following. The end mirrors reflect the beams 
such that they are recombincd at the beam splitter. The South and West ports of the beam splitter 
are possible output port, however, in many cases, only the South port is used. 



5.2 Michelson interferometer 



We came across the Michelson interferometer in Section 2.4 when we discussed the phase relation 
at a beam splitter. The typical optical layout of the Michelson interferometer is shown again in 
Figure |28j a laser beam is split by a beam splitter and send along two perpendicular interferometer 
arms. The four directions seen from the beam splitter are called North, East, West and South. The 
ends of these arms (North and East) are marked by highly reflective end mirrors, which reflect the 
beams back into themselves so that they can be recombined by the beam splitter. Generally, the 
Michelson interferometer has two outputs, namely the so far unused beam splitter port (South) and 
the input port (West). Both output ports can be used to obtain interferometer signals, however, 
most setups are designed such that the signals with high signal-to-noise ratios are detected in the 
South port. 

The Michelson interferometer output is determined by the laser wavelength A, the reflectivity 
and transmittance of the beam splitter and the end mirrors, and the length of the interferometer 
arms. In many cases the end mirrors are highly reflective and the beam splitter ideally a 50:50 
beam splitter. In that case, we can compute the output for a monochromatic field as shown in 



Section 2.4 Using Equation (20 1 we can write the field in the South port as 



Es = Eo'-{e'^''^'' + e'2^^^). (96) 



We define the common arm length and the arm-length difference as 

L = ' 



which yield 2Ljv = 2L + AL and 2Le = 2L — AL. Thus, we can further simplify to get 

Es = Eo '-e'^^^ (e"=^^ + e"'^^^) = Eo ie'^'^-^ cos(fcAL). (98) 
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The photo detector then produces a signal proportional to 



S = EsE*s = Pq cos^ikAL) = Pq cos^ {2tt AL / X) . 



(99) 



This signal is depicted in Figure [29j it shows that the power in the South port changes between 
zero and the input power with a period of AL/A = 0.5. The tuning at which the output power 
drops to zero is called the dark fringe. Current interferometric gravitational- wave detectors operate 
their Michelson interferometer at or near the dark fringe. 




The above seems to indicate that the macroscopic arm-length difference plays no role in the 
Michelson output signal. However, this is only correct for a monochromatic laser beam with infinite 
coherence length. In real interferometers care must be taken that the arm-length difference is well 
below the coherence length of the light source. In gravitational-wave detectors the macroscopic 
arm-length difference is an important design feature; it is kept very small in order to reduce 
coupling of laser noise into the output but needs to retain a finite size to allow the transfer of 
phase modulation sidebands from the input to the output port; this is illustrated in the Finesse 
example below and will be covered in detail in Section [6. 4[ 
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5.3 Finesse examples 
5.3.1 Michelson power 



fexample_michelson_power Wed Aug 19 15:01 :57 2009 




50 100 150 200 250 300 

north mirror tuning [deg] 



Figure 30: Finesse example: Michelson power. 

The power in the South port of a Michelson detector varies as the cosine squared of the microscopic 
arm length difference. The maximum output can be equal to the input power, but only if the 
Michelson interferometer is symmetric and lossless. The tuning for which the South port power is 
zero is referred to as the dark fringe. 

Finesse input file for 'Michelson power' 

laser 11 1 nl '/, laser with P=1W at the default frequency 

space si 1 1 nl n2 space of Im length 

runl: bs bl 0.5 0.5 n2 nNl nEl nSl 7, 50:50 beam splitter 
run2: bs bl 0.4 0.6 n2 nNl nEl nSl 40:60 beam splitter 
runS: bs bl 0.45 0.45 n2 nNl nEl nSl 7. 45:45 beam splitter 
space LN 1 1 nNl nN2 7, north arm 
space LE 1 1 nEl nE2 7> east arm 

mirror mN 1 nN2 dump 7, north end mirror, lossless 
mirror mE 1 nE2 dump 7. east end mirror, lossless 
space s2 1 1 nSl nout 

pd South nout 7 photo detector in South port 

xaxis mN phi lin 300 100 7 changing the microscopic position of mN 



5.3.2 Michelson modulation 

This example demonstrates how a macroscopic arm length difference can cause different 'dark 
fringe' tuning for injected fields with different frequencies. In this case, some of the 10 MHz 
modulation sidebands are transmitted when the interferometer is tuned to a dark fringe for the 
carrier light. This effect can be used to separate light fields of different frequencies. It is also the 
cause for transmission of laser noise (especially frequency noise) into the Michelson output port 
when the interferometer is not perfectly symmetric. 
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fexample_michelson_modulation Wed Aug 19 16:49:14 2009 
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Figure 31: Finesse example: Michelson modulation. 



Finesse input file for 'Michelson modulation' 

laser 11 1 nl laser with P=1W at the default frequency 

space si 1 1 nl n2 '/, space of Im length 

mod eoml lOM 0.3 1 pm n2 n3 phase modulation at 10 MHz 

space s2 1 1 nS n4 another space of Im length 

bs bl 0.5 0.5 n4 nNl nEl nSl 50:50 beam splitter 

space LN 100 1 nNl nN2 north arm 

space LE 110 1 nEl nE2 7, east arm 

mirror mN 1 nN2 dump 7, north end mirror, lossless 
mirror mE 1 nE2 dump 7 east end mirror, lossless 
space s3 1 1 nSl nout 

ad carrier nout 7 amplitude detector for carrier field 

ad sideband lOM nout 7 amplitude detector for +10 MHz sideband 

xaxis mN phi lin 300 100 7 changing the microscopic position of mN 
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6 Interferometric Length Sensing and Control 



In this section we introduce interferometers as length sensing devices. In particular, we explain how 
the Fabry-Perot interferometer and the Michelson interferometer can be used for high-precision 
measurements and that both require a careful control of the base length (which is to be measured) 
in order to yield their large sensitivity. In addition, we briefly introduce the general concepts of 
error signals and transfer functions, which are used to describe most essential features of length 
sensing and control. 



6.1 Error signals and transfer functions 

In general, we will call an error signal any measured signal suitable for stabilising a certain ex- 
perimental parameter p with a servo loop. The aim is to maintain the variable p at a user-defined 
value, the operating point, pq. Therefore, the error signal must be a function of the parameter p. 
In most cases it is preferable to have a bipolar signal with a zero crossing at the operating point. 
The slope of the error signal at the operating point is a measure of the 'gain' of the sensor (which 
in the case of interferometers is a combination of optics and electronics) . 

Transfer functions describe the propagation of a periodic signal through a plant and are usually 
given as plots of amplitude and phase over frequency. By definition a transfer function describes 
only the linear coupling of signals inside a system. This means a transfer function is independent 
of the actual signal size. For small signals or small deviations, most systems can be linearised and 
correctly described by transfer functions. 

Experimentally, network analysers are commonly used to measure a transfer function: one 
connects a periodic signal (the source) to an actuator of the plant (which is to be analysed) and to 
an input of the analyser. A signal from a sensor that monitors a certain parameter of the plant is 
connected to the second analyser input. By mixing the source with the sensor signal the analyser 
can determine the amplitude and phase of the input signal with respect to the source (amplitude 
equals one and the phase equals zero when both signals are identical). 

Mathematically, transfer functions can be modeled similarly: applying a sinusoidal signal 
sin(aj<,t) to the interferometer, e.g., as a position modulation of a cavity mirror, will create phase 
modulation sidebands with a frequency offset of ±cjs to the carrier light. If such light is detected 
in the right way by a photodiode, it will include a signal at the frequency component ujs, which 



can be extracted, for example, by means of demodulation (see Section 4.2). 

Transfer functions are of particular interest in relation to error signals. Typically a transfer 
function of the error signal is required for the design of the respective electronic servo. A 'transfer 
function of the error signal' usually refers to a very specific setup: the system is held at its operating 
point, such that, on average, p — pq. A signal is applied to the system in the form of a very small 
sinusoidal disturbance of p. The transfer function is then constructed by computing for each signal 
frequency the ratio of the error signal and the injected signal. Figure [32] shows an example of an 
error signal and its corresponding transfer function. The operating point shall be at 

Xd^O and a;Ep(xd = 0)=0 (100) 

The optical transfer function Topt,Xd with respect to this error signal is defined by 

^Ep(/) = 7'opt,Xd'7dot^d(/), (101) 

with Tdet as the transfer function of the sensor. In the following, Tjot is assumed to be unity. At 
the zero crossing the slope of the error signal represents the magnitude of the transfer function for 
low frequencies: 

'^'^^^ , =\To,u^.\y^, (102) 
Ud=o 



dxd 
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Figure 32: Example of an error signal: the top graph shows the electronic interferometer output 
signal as a function of mirror displacement. The operating point is given as the zero crossing, and 
the error-signal slope is defined as the slope at the operating point. The right graph shows the 
magnitude of the transfer function mirror displacement — error signal. The slope of the error 
signal (left graph) is equal to the low frequency limit of the transfer function magnitude (see 
Equation (Il02])). 
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The quantity above will be called the error-signal slope in the following text. It is proportional to 
the optical gain iTopt,^;^!, which describes the amplification of the gravitational- wave signal by the 
optical instrument. 
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Figure 33: The top plot shows the cavity power as a function of the cavity tuning. A tuning of 
360° refers to a change in the cavity length by one laser wavelength. The bottom plot shows the 
differentiation of the upper trace. This illustrates that near resonance the cavity power changes 
very rapidly when the cavity length changes. However, for most tunings the cavity seems not 
sensitive at all. 



6.2 Fabry— Perot length sensing 

In Figure |25] we have plotted the circulating power in a Fabry-Perot cavity as a function of the 
laser frequency. The steep features in this plot indicate that such a cavity can be used to measure 



changes in the laser frequency. From the equation for the circulating power (see Equation (88)) 



= l + i?ii?2-2rir2COs(2/fcL) = ^^"^^ 

we can see that the actual frequency dependence is given by the cos(2fcL) term. Writing this term 
as 

cos{2kL) = cos (^27r— ^ , (104) 

we can highlight the fact that the cavity is in fact a reference for the laser frequency in relation to 
the cavity length. If we know the cavity length very well, a cavity should be a good instrument to 
measure the frequency of a laser beam. However, if we know the laser frequency very accurately, 
we can use an optical cavity to measure a length. In the following we will detail the optical setup 
and behaviour of a cavity used for a length measurement. The same reasoning applies for frequency 
measurements. If we make use of the resonant power enhancement of the cavity to measure the 
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Figure 34: Typical setup for using the Pound-Drever-Hall scheme for length sensing and with a 
two- mirror cavity: the laser beam is phase modulated with an electro-optical modulator (EOM). 
The modulation frequency is often in the radio frequency range. The photodiode signal in reflection 
is then electrically demodulated at the same frequency. 



cavity length, we can derive the sensitivity of the cavity from the differentiation of Equation 
which gives the slope of the trace shown in Figure [25j 

dPi/Po -4Tirir2fcsin(2fcL) 

= ' 



with d as defined in Equation ( 103 1 . This is plotted in Figure 33 together with the cavity power 



as a function of the cavity tuning. From Figure [33] we can deduce a few key features of the cavity: 

• The cavity must be held as near as possible to the resonance for maximum sensitivity. This 
is the reason that active servo control systems play an important role in modern laser inter- 
ferometers. 

• If we want to use the power directly as an error signal for the length, we cannot use the 
cavity directly on resonance because there the optical gain is zero. A suitable error signal 
(i.e., a bipolar signal) can be constructed by adding an offset to the light power signal. A 
control system utilising this method is often called DC-lock or offset-lock. However, we show 
below that more elegant alternative methods for generating error signals exist. 

• The differentiation of the cavity power looks like a perfect error signal for holding the cavity on 
resonance. A signal proportional to such differentiation can be achieved with a modulation- 
demodulation technique. 



6.3 The Pound— Drever— Hall length sensing scheme 

This scheme for stabilising the frequency of a light field to the length of a cavity, or vice versa, is 
based on much older techniques for performing very similar actions with microwaves and microwave 
resonators. Drever and Hall have adapted such techniques for use in the optical regime [IS] and 
today what is now called the Pound-Drever-Hall technique can be found in a great number of 
different types of optical setups. An example layout of this scheme is shown in Figure [34j in this 
case for generating a length (or frequency) signal of a two-mirror cavity. The laser is passed through 
an electro-optical modulator, which applies a periodic phase modulation at a fixed frequency. In 
many cases the modulation frequency is chosen such that it resides in the radio frequency band for 
which low-cost, low-noise electronic components are available. The phase modulated light is then 
injected into the cavity. However, from the frequency domain analysis introduced in Section [5] 
we know that in most cases not all the light can be injected into the cavity. Let's consider the 
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Figure 35: This figure shows an example of a Pound-Drever-HaU (PDH) signal of a two-mirror 
cavity. The plots refer to a setup in which the cavity mirrors are stationary and the frequency of 
the input laser is tuned linearly. The upper trace shows the light power circulating in the cavity. 
The three peaks correspond to the frequency tunings for which the carrier (main central peak) or 
the modulation sidebands (smaller side peaks) are resonant in the cavity. The lower trace shows 
the PDH signal for the same frequency tuning. Coincident with the peaks in the upper trace 
are bipolar structures in the lower trace. Each of the bipolar structures would be suitable as a 
length-sensing signal. In most cases the central structure is used, as experimentally it can be easily 
identified because its slope has a different sign compared to the sideband structures. 
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example of an over-coupled cavity with the reflectivity of the end mirror R2 < 1. Such a cavity 



would have a frequency response as shown in the top traces of Figure 26 (recall that the origin 
of the frequency axis refers to an arbitrarily chosen default frequency, which for this figure has 
been selected to be a resonance frequency of the cavity) . If the cavity is held on resonance for the 
unmodulated carrier field, this field enters the cavity, gets resonantly enhanced and a substantial 
fraction is transmitted. If the frequency offset of the modulation sidebands is chosen such that it 
does not coincide with (or is near to) an integer multiple of the cavity's free spectral range, the 
modulation sidebands are mostly reflected by the cavity and will not be influenced as much by 
the resonance condition of the cavity as the carrier. The photodiode measuring the reflected light 
will see the optical beat between the carrier fleld and the modulation sidebands. This includes a 
component at the modulation frequency which is a measure of the phase difference between the 
carrier field and the sidebands (given the setup as described above). Any slight change of the 
cavity length would introduce a proportional change in the phase of the carrier field and no change 
in the sideband fields. Thus the photodiode signal can be used to measure the length changes of 
the cavity. One of the advantages of this method is the fact that the so-generated signal is bipolar 
with a zero crossing and steep slope exactly at the cavity's resonance, see Figure [35} 
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Figure 36: Power and slope of a Michelson interferometer. The upper plot shows the output power 



of a Michelson interferometer as detected in the South port (as already shown in Figure 29 1 . The 



lower plot shows the optical gain of the instrument as given by the slope of the upper plot. 



6.4 Michelson length sensing 

Similarly to the two-mirror cavity, we can start to understand the length-sensing capabilities of the 
Michelson interferometer by looking at the output light power as a function of a mirror movement, 
as shown in Figure [29] The power changes as sine squared with the maximum slope at the point 
when the output power (in what we call the South port) is half the input power. The slope of the 
output power, which is the optical gain of the instrument for detecting a differential arm-length 



48 



change AL with a photo detector in the South port can be written as 



dAL A V A 



(106) 



and is shown in Figure [36| The most notable difference of the optical gain of the Michelson 



interferometer with respect to the Fabry-Perot interferometer (see Figure 33) is the wider, more 
smooth distribution of the gain. This is due to the fact that the cavity example is based on a high- 
finesse cavity in which the optical resonance effect is dominant. In a basic Michelson interferometer 
such resonance enhancement is not present. 

However, the main difference is that the measurement is made differentially by comparing 
two lengths. This allows one to separate a larger number of possible noise contributions, for 
example noise in the laser light source, such as amplitude or frequency noise. This is why the 
main instrument for gravitational-wave measurements is a Michelson interferometer. However, the 
resonant enhancement of light power can be added to the Michelson, for example, by using Fabry- 
Perot cavities within the Michelson. This construction of new topologies by combining Michelson 
and Fabry-Perot interferometers will be described in detail in a future version of this review. 

The Michelson interferometer has two longitudinal degrees of freedom. These can be represented 
by the positions (along the optical axes) of the end mirrors. However, it is more efhcient to use 
proper linear combinations of these and describe the Michelson interferometer length or position 



information by the common and differential arm length, as introduced in Equation ( 97 1 : 



AL = Ln - L 



The Michelson interferometer is intrinsically insensitive to the common arm length L. 
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Figure 37: This length sensing scheme is often referred to as frontal or Schnupp modulation: an 
EOM is used to phase modulate the laser beam before entering the Michelson interferometer. The 
signal of the photodiode in the South port is then demodulated at the same frequency used for the 
modulation. 



6.5 The Schnupp modulation scheme 

Similar to the Fabry-Perot cavity, the Michelson interferometer is also often used to set an oper- 
ating point where the optical gain of a direct light power detection is zero. This operating point, 
given by AL/X = {2N + 1) • 0.25 with N a non-negative integer, is called dark fringe. This oper- 
ating point has several advantages, the most important being the low (ideally zero) light power on 
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the diode. Highly efficient and low-noise photodiodes usually use a small detector area and thus 
are typically not able to detect large power levels. By using the dark fringe operating point, the 
Michelson interferometer can be used as a null instrument or null measurement, which generally 
is a good method to reduce systematic errors [35j. 

One approach to make use of the advantages of the dark fringe operating point is to use an 
operating point very close to the dark fringe at which the optical gain is not yet zero. In such a 
scenario a careful trade-off calculation can be done by computing the signal-to-noise with noises 
that must be suppressed, such as the laser amplitude noise. This type of operation is usually 
referred to as DC control or offset control and is very similar to the similarly-named mechanism 
used with Fabry-Perot cavities. 

Another option is to employ phase modulated light, similar to the Pound-Drever-Hall scheme 
described in Section |6.3[ The optical layout of such a scheme is depicted in Figure [37} an electro- 
optical modulator is used to apply a phase modulation at a fixed (usually RF type) frequency to 
the (monochromatic) laser light before it enters the interferometer. The photodiode signal from 
the interferometer output is then demodulated at the same frequency. This scheme allows one to 
operate the interferometer precisely on the dark fringe. The method originally proposed by Lise 
Schnupp is also sometimes referred to as frontal modulation. 

The optical gain of a Michelson interferometer with Schnupp modulation is shown in Figure [39] 
in Section [631 
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6.6 Finesse examples 
6.6.1 Cavity power and slope 
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Figure 38: Finesse example: Cavity power and slope. 

Figure |33] shows a plot of the analytical functions describing the power inside a cavity and its 
differentiation by the cavity tuning. This example recreates the plot using a numerical model in 
Finesse. 



Finesse input file for 'Cavity power and slope' 



laser 11 1 nl 7, laser with P=1W at the default frequency 

space si 1 1 nl n2 7, space of Im length 

mirror ml . 9 . 1 n2 n3 7. cavity input mirror 

space L 1200 1 nS n4 7. cavity length of 1200m 

mirror m2 1.0 0.0 n4 dump 7« cavity output mirror 

pd P n3 7. photo diode measuring the intra-cavity power 

7, for the plot we perform two sequenctial runs of Finesse using 'mkat' 
7« 1) first trace: plot the power (switching to log plot) 
runl : yaxis log abs 

7. 2) second trace; plot the differentiation 
run2: diff m2 phi 

xaxis m2 phi lin -50 250 300 7o changing the microscopic tuning of mirror m2 
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6.6.2 Michelson with Schnupp modulation 
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Figure 39: Finesse example: Michelson with Schnupp modulation. 

Figure |39] shows the demodulated photodiode signal of a Michelson interferometer with Schnupp 
modulation, as well as its differentiation, the latter being the optical gain of the system. Comparing 
this figure to Figure |36[ it can be seen that with Schnupp modulation, the optical gain at the dark 
fringe operating points is maximised and a suitable error signal for these points is obtained. 

Finesse input file for 'Michelson with Schnupp modulation' 



laser 11 1 nl laser with P=1W at the default frequency 

space si 1 1 nl n2 7, space of Im length 

mod eoml lOM 0.3 1 pm n2 n3 7. phase modulation at 10 MHz 

space s2 1 1 nS n4 7, another space of Im length 

bs bl 0.5 0.5 n4 nNl nEl nSl 7. 50:50 beam splitter 

space LN 100 1 nNl nN2 % north arm 

space LE 110 1 nEl nE2 7 east arm 

mirror mN 1 22 nN2 dump 7 north end mirror, lossless 
mirror mE 1 -22 nE2 dump 7 east end mirror, lossless 
space s3 1 1 nSl nout 



runl : pdl South lOM -115 nout 7 demodulated output signal 
run2: pdl South lOM -115 nout 7 demodulated output signal 
run2: diff mN phi 7 computing the slope of the signal 

xaxis mN phi lin 300 100 7 changing the microscopic position of mN 

put mE phi $mxl 7 moving mE as -mN to make a differential motion 
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7 Beam Shapes: Beyond the Plane Wave Approximation 



In previous sections we have introduced a notation for describing the on-axis properties of electric 
fields. Specifically, we have described the electric fields along an optical axis as functions of 
frequency (or time) and the location z. Models of optical systems may often use this approach for 
a basic analysis even though the respective experiments will always include fields with distinct off- 
axis beam shapes. A more detailed description of such optical systems needs to take the geometrical 
shape of the light field into account. One method of treating the transverse beam geometry is to 
describe the spatial properties as a sum of 'spatial components' or 'spatial modes' so that the 
electric field can be written as a sum of the different frequency components and of the different 
spatial modes. Of course, the concept of modes is directly related to the use of a sort of oscillator, 
in this case the optical cavity. Most of the work presented here is based on the research on laser 
resonators reviewed originally by Kogelnik and Li [35^ . Siegman has written a very interesting 
historic review of the development of Gaussian optics [5T1 [SD] and we use whenever possible the 
same notation as used in his textbook 'Lasers' [52]- 

This section introduces the use of Gaussian modes for describing the spatial properties along 
the transverse orthogonal x and y directions of an optical beam. We can write 

E{t,x,y,z) ^ a-jnm Unmix, y,z) ex.p{i{u!jt - kjz)), (107) 

j n.m 

with Unm ELS special functions describing the spatial properties of the beam and aj^m as complex 
amplitude factors (wj is again the angular frequency and kj — ujj/c). For simplicity we restrict 
the following description to a single frequency component at one moment in time (t = 0), so 

E{x,y,z) ^ exp{~ikz) ^ anm Unm{x,y, z). (108) 

In general, different types of spatial modes Unm can be used in this context. Of particular interest 
are the Gaussian modes, which will be used throughout this document. Many lasers emit light that 
closely resembles a Gaussian beam: the light mainly propagates along one axis, is well collimated 
around that axis and the cross section of the intensity perpendicular to the optical axis shows a 
Gaussian distribution. The following sections provide the basic mathematical framework for using 
Gaussian modes for analysing optical systems. 



7.1 The paraxial wave equation 

Mathematically, Gaussian modes are solutions to the paraxial wave equation - a specific wave 
equation for electromagnetic fields. All electromagnetic waves are solutions to the general wave 
equation, which in vacuum can be given as: 

A£;-^£; = G. (109) 

But laser light fields are special types of electromagnetic waves. For example, they are characterised 
by low diffraction. Hence, a laser beam will have a characteristic length w describing the 'width' 
(the dimension of the field transverse to the main propagation axis) , and a characteristic length I 
defining some local length along the propagation over which the beam characteristics do not vary 
much. By definition, for what we call a beam w is typically small and I large in comparison, so 
that w/l can be considered small. In fact, the paraxial wave equation (and its solutions) can be 



derived as the first-order terms of a series expansion of Equation (109) into orders of w/l j38j 
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A simpler approach to the paraxial-wave equation goes as follows: A particular beam shape 
shall be described by a function u{x, y, z) so that we can write the electric field as 

E{x,y, z) = u(x,y, z) exp(— ifcz). (HO) 

Substituting this into the standard wave equation yields a differential equation for u: 

{dl + dl+ dl) u{x, y, z) - 2i fc9^u(x, y, z) = 0. (Ill) 

Now we put the fact that u(x^ y, z) should be slowly varying with z in mathematical terms. The 
variation of u{x, y, z) with z should be small compared to its variation with x or y. Also the second 
partial derivative in z should be small. This can be expressed as 

\dlu{x,y,z)\ < \2kdMx,y,z)\ , \dlu{x,y,z)\ , \d'^u{x,y, z)\ . (112) 



With this approximation, Equation (111 I can be simplified to the paraxial wave equation, 

{dl + dl) u{x, y, z) - 2i kdM^, y, z) = 0. (113) 
Any field u that solves this equation represents a paraxial beam shape when used in the form given 



in Equation ( 1 10 1 



7.2 Transverse electromagnetic modes 



In general, any solution u{x,y,z) of the paraxial wave equation. Equation (113), can be employed 
to represent the transverse properties of a scalar electric field representing a beam-like electro- 
magnetic wave. Especially useful in this respect are special families or sets of functions that are 
solutions of the paraxial wave equation. When such a set of functions is complete and countable, 
it's called a set of transverse electromagnetic modes (TEM). For instance, the set of Hermite- 
Gauss modes are exact solutions of the paraxial wave equation. These modes are represented 
by an infinite, countable and complete set of functions. The term complete means they can be 
understood as a base system of the function space defined by all solutions of the paraxial wave 
equation. In other words, we can describe any solution of the paraxial wave equation u' by a linear 
superposition of Hermite-Gauss modes: 

u'{x,y,z) = a J 

n.m 

which in turn allows us to describe any laser beam using a sum of these modes: 

E{t,x,y,z)^Y Y a-jnm Unm{x,y,z) expiiiujjt - kjZ)). (115) 



The Hermite-Gauss modes as given in this document (see Section 7.5) are orthonormal so that 



* flifn^n' and m = m' 1 , , 

dxdyu„„«„,,„,=.<5„„,5„„, = | P otherwise J' ^^^^^ 

This means that, in the function space defined by the paraxial wave equation, the Hermite-Gauss 
functions can be understood as a complete set of unit-length basis vectors. This fact can be 
utilised for the computation of coupling factors. Furthermore, the power of a beam, as given by 



Equation (108), being detected on a single-element photodetector (provided that the area of the 



detector is large with respect to the beam) can be computed as 

EE* ^Y.°'^^^<rn, (117) 
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or for a beam with several frequency components (compare with Equation (76l) as 

= mm a™ma*„™ with I J e {0, . . . , TV} A u, = u,}. (118) 
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Figure 40: One dimensional cross-section of a Gaussian beam. The width of the beam is given by 
the radius w at which the intensity is 1 /e^ of the maximum intensity. 



7.3 Properties of Gaussian beams 

The basic or 'lowest-order' Hermite-Gauss mode is equivalent to what is usually called a Gaussian 
beam and is given by 



1 



TT w{z) 



exp(i5'(z)) exp — ifc 



y 



y 



2Rc{z) w-^{z) 



(119) 



The parameters of this equation are explained in detail below. The shape of a Gaussian beam is 
quite simple: the beam has a circular cross section, and the radial intensity profile of a beam with 
total power P is given by 

op 

^(^) = (-2rVii;2) , (120) 

with w the spot size, defined as the radius at which the intensity is 1/e^ times the maximum 



intensity /(O). This is a Gaussian distribution, see Figure 40 hence the name Gaussian beam. 

Figure [41] shows a different cross section through a Gaussian beam: it plots the beam size as a 
function of the position on the optical axis. 

Such a beam profile (for a beam with a given wavelength A) can be completely determined by 
two parameters: the size of the minimum spot size wq (called beam waist) and the position zq of 
the beam waist along the 2-axis. 

To characterise a Gaussian beam, some useful parameters can be derived from wq and zq. A 
Gaussian beam can be divided into two different sections along the z-axis: a near field - a region 
around the beam waist, and a far field - far away from the waist. The length of the near-field 
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near field ' far field 

Figure 41: Gaussian beam profile along z: this cross section along the a;-z-plane illustrates how 
the beam size w{z) of the Gaussian beam changes along the optical axis. The position of minimum 
beam size wq is called beam waist. See text for a description of the parameters 0, zji and Re- 



region is approximately given by the Rayleigh range zr. The Rayleigh range and the spot size are 
related by 

zn-"^. (121) 
With the Rayleigh range and the location of the beam waist, we can usefully write 

w{z) = wo]jl+(^^^^. (122) 

This equation gives the size of the beam along the z-axis. In the far-field regime (z ^ zr, Zq), it 
can be approximated by a linear equation, when 

/ \ z z\ . 
u;(z wuJo— = . (123) 

The angle Q between the z-axis and •w{z) in the far field is called the diffraction angl^anA is 
defined by 

Q — arctan ( — ) = arctan ( ) sa — . (124) 

\ZkJ \ TTtWo / 2r 

Another useful parameter is the radius of curvature of the wavefront at a given point z. The 
radius of curvature describes the curvature of the 'phase front' of the electromagnetic wave - a 
surface across the beam with equal phase - intersecting the optical axis at the position z. We 

^Also known as the far-field angle or the divergence of the beam. 
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obtain the radius of curvature as a function of z: 



We also find: 



Rc{z)^z-zo + ^^. (125) 

Z- Zq 

Rc « oo, z — zo ^ (beam waist) 

Rc ^ z, z > zr, zo (far field) (126) 

Rc — 2zr, z — zq = zr (maximum curvature). 

7.4 Astigmatic beams: the tangential and sagittal plane 

If the interferometer is confined to a plane (here the x-z plane), it is convenient to use projections 
of the three-dimensional description into two planes [IBj: the tangential plane, defined as the x-z 
plane and the sagittal plane as given by y and z. 

The beam parameters can then be split into two respective parameters: zg^s, wo.s for the sagittal 
plane and zq,* and wo,* for the tangential plane so that the Hcrmite-Gauss modes can be written 
as 

Unm{x, V) = U„{X, Zo,t, Wo,t) Umiu, Zq^s, Wo,s)- (127) 

Beams with different beam waist parameters for the sagittal and tangential plane are astigmatic. 

Remember that these Hermite-Gauss modes form a base system. This means one can use 
the separation into sagittal and tangential planes even if the actual optical system does not show 
this special type of symmetry. This separation is very useful in simplifying the mathematics. In 
the following, the term beam parameter generally refers to a simple case where wo,x — WQ,y and 
zo,x = ZQ^y but all the results can also be applied directly to a pair of parameters. 

7.5 Higher-order Hermite— Gauss modes 

The complete set of Hcrmite-Gauss modes is given by an infinite discrete set of modes Un,„(a;, y, z) 
with the indices n and m as mode numbers. The sum n+m is called the order of the mode. The 
term higher-order modes usually refers to modes with an order n + m > 0. The general expression 
for Hermite-Gauss modes can be given as [55] 

Mnm(a;, y, Z) = Unix, z)Um{y, z) , (128) 

with 



, (t A- f 2 '1 1/4 / cxp(i(2n+l)M^(^)) \^/^ 
ln[-l, Z) — y^) y 2"n!u)(z) J 



(129) 



and Hn{x) the Hermite polynomials of order n. The first Hermite polynomials, without normali- 
sation, can be written 

Ho{x) = l Hi{x)^2x 
H2{x)= 4x^-2 Halx) ^ Sx^ - I2x. 

Further orders can be computed recursively since 



(130) 



Hn+i{x) ^ 2xHnix) - 2nH„^i{x). (131) 
For both transverse directions we can also rewrite the above to 

-1/2 



Unmix, y,z) = (2"+™~in!m!7r) ' exp (i (n + m + l)^'(z)) x 



(132) 
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The latter form has the advantage of clearly showing the extra phase shift along the z-axis of 
{n + m + l)5'(z), called the Gouy phase; see Section 7.8 




m 1 




Figure 42: This plot shows the intensity distribution of Hermite-Gauss modes One can 

see that the intensity distribution becomes wider for larger mode indices and the peak intensity 
decreases. The mode index defines the number of dark stripes in the respective direction. 



7.6 The Gaussian beam parameter 

For a more compact description of the interaction of Gaussian modes with optical components we 
will make use of the Gaussian beam parameter q 36 . The beam parameter is a complex quantity 
defined as 

' ' ■ ' (133) 



It can also be written as 



q{z) = i zr, + z — zo = go + 2 — and go = i z-r- (134) 



Using this parameter, Equation (119) can be rewritten as 

u{x,y,z) = -^ex]^(-\k ^^\ V (135) 
g(z) V 2g(z) J 

Other parameters, like the beam size and radius of curvature, can also be written in terms of the 
beam parameter q: 

.^W = ^ifl. (136) 
»S = (137) 
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HGOO amplitude 



HGOO phase 






Figure 43: These plots show the amphtude distribution and wave front (phase distribution) of 
Hermite-Gaussian modes Unm (labeled as HGnm in the plot). All plots refer to a beam with 
A = 1 iim, w = 1 mm and distance to waist z = 1 m. The mode index (in one direction) defines 
the number of zero crossings (along that axis) in the amplitude distribution. One can also see that 
the phase distribution is the same spherical distribution, regardless of the mode indices. 
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ZR^^{q} (138) 

and 

«c(.) ^ (139) 

The Hermite-Gauss modes can also be written using the Gaussian beam parameter a^ 

Unmix, y, z) = Unix, z)u„,iy , z) with 

u(xz)- i^A^'^ i^A r *^r^' H exp f-i ^^^^^ 

7.7 Properties of higher-order Hermite— Gauss modes 

Some of the properties of Hermite-Gauss modes can easily be described using cross sections of the 
field intensity or field amplitude. Figure [42| shows such cross sections, i.e., the intensity in the x-y 
plane, for a number of higher-order modes. This shows a x-y symmetry for mode indices n and 
m. We can also see how the size of the intensity distribution increases with the mode index, while 
the peak intensity decreases. 



Similarly, Figure 43 shows the amplitude and phase distribution of several higher-order Hermite- 



Gauss modes. Some further features of Hermite-Gauss modes: 

• The size of the intensity profile of any sum of Hermite-Gauss modes depends on z while its 
shape remains constant over propagation along the optical axis. 

• The phase distribution of Hermite-Gauss modes shows the curvature (or radius of curvature) 
of the beam. The curvature depends on z but is equal for all higher-order modes. 

Note that these are special features of Gaussian beams and not generally true for arbitrary 
beam shapes. Figure [44| for example, shows the amplitude and phase distribution of a triangular 
beam at the point where it is (mathematically) created and after a 10 m propagation. Neither the 
shape is preserved nor does it show a spherical phase distribution. 

7.8 Gouy phase 

The equation for Hermite-Gauss modes shows an extra longitudinal phase lag. This Gouy phase 
[27l[26] describes the fact that, compared to a plane wave, the Hermite-Gauss modes have a slightly 
slower phase velocity, especially close to the waist. The Gouy phase can be written as 

*(z) arctan ( \ ^ (141) 



ZR 



or, using the Gaussian beam parameter. 



l'(z) = arctan I §^ ) . (142) 



Compared to a plane wave, the phase lag (/? of a Hermite-Gauss mode is 

(n + m + l)*(z). (143) 



^Please note that this formula from 52 is very compact. Since the parameter g is a complex number, the 
expression contains at least two complex square roots. The complex square root requires a different algebra than 
the standard square root for real numbers. Especially the third and fourth factors can not be simplified in any 

obvious wav ( ™ ( qoq*(^) Y'^ ^ ( 1o+^l'^(- )V^ 

Obvious way. V^^J^j ^ U"+i(z),*' 



I 
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Amplitude Phase 




Amplitude Phase 




Figure 44: These top plots show a triangular beam shape and phase distribution and the bottom 
plots the diffraction pattern of this beam after a propagation of 2 = 5 m. It can be seen that the 
shape of the triangular beam is not conserved and that the phase front is not spherical. 
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With an astigmatic beam, i.e., different beam parameters in the tangential and sagittal planes, 
this becomes 

f=(n+-j-i>t{z)+(m+-j^s{z), (144) 



with 

as the Gouy phase in the tangential plane (and ^'^ is similarly defined in the sagittal plane). 



*,(z) = arctan(|M ) , (I45) 



7.9 Laguerre— Gauss modes 

Laguerre-Gauss modes are another complete set of functions, which solve the paraxial wave equa- 
tion. They are defined in cylindrical coordinates and can have advantages over Hermite-Gauss 
modes in the presence of cylindrical symmetry. More recently, Laguerre-Gauss modes are being 
investigated in a different context: using a pure higher-order Laguerre-Gauss mode instead of the 
fundamental Gaussian beam can significantly reduce the impact of mirror thermal noise on the 
sensitivity of gravitational wave detectors [SH [T2] . Laguerre-Gauss modes are commonly given 
as [32] 

^pAr,^,z) = 4i)^^^feexp(i(2p-H |Z| -I- l)*(z)) 

X (^)^r(4S-)exp(-ifc-^+i^0), 

with r, (j) and z as the cylindrical coordinates around the optical axis. The letter p is the radial 
mode index, I the azimuthal mode indejj^ and Lp ^ (x) are the associated Laguerre polynomials: 

w 4 |:J(;t^ )<-)'■ a") 

All other parameters {■w{z),q{z), . . .) are defined as above for the Hermite-Gauss modes. 



The dependence of the Laguerre modes on <j) as given in Equation ( 146 1 results in a spiraling 
phase front, while the intensity pattern will always show unbroken concentric rings; see Figure [45] 
These modes are also called helical Laguerre-Gauss modes because of the their special phase 
structure. 

The reader might be more familiar with a slightly different type of Laguerre modes (compare 



Figure 46 and Figure 47 1 that features dark radial lines as well as dark concentric rings. Mathe- 



matically, these can be described simply by replacing the phase factor exp(i/(/)) in Equation (146) 
by a sine or cosine function. For example, an alternative set of Laguerre-Gauss modes is given 
by [SS] 

K%r,c^,z) = ^V' (WI^I+P)I ^>^P(U2P+I^I + W)) 

This type of mode has a spherical phase front, just as the Hermite-Gauss modes. We will refer to 
this set as sinusoidal Laguerre-Gauss modes throughout this document. 

For the purposes of simulation it can be sometimes useful to decompose Laguerre-Gauss modes 
into Hermite-Gauss modes. The mathematical conversion for helical modes is given as [71 [T] 

N 

u'^%{x,y,z) = Y,'^%{n,m,k)u%%^i^{x,y,z), (149) 

fc=0 



|52| states that the indices must obey the following relations: < |/| < p. However, that is not the case. 
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Figure 45: These plots show the amphtude distribution and wave front (phase distribution) of 
hehcal Laguerre-Gauss modes Upi. All plots refer to a beam with A = 1 pm, w = 1 mm and 
distance to waist 2 = 1 m. 
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with real coefficients 

if = n + m. This relates to the common definition of Laguerre modes as Upi as follows: p — 
min(n, m) and I = n — m. The coefficients to, fc) can be computed numerically by using Jacobi 
polynomials. Jacobi polynomials can be written in various forms: 

Pn^i^) = - + ^r^id.ni - + ^f'"', (151) 

or 

K'^i-) = ^ E )(- - + (152) 

which leads to 

= \/^^Sr (-2)'=^^r''"-'(o). (153) 

I — 
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Figure 46: Intensity profiles for helical Laguerre-Gauss modes Upi. The uqq mode is identical to 
the Hermite-Gauss mode of order 0. Higher-order modes show a widening of the intensity and 
decreasing peak intensity. The number of concentric dark rings is given by the radial mode index 
P- 



7.10 Tracing a Gaussian beam through an optical system 

Whenever Gauss modes are used to analyse an optical system, the Gaussian beam parameters (or 
equivalent waist sizes and locations) must be defined for each location at which field amplitudes are 
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Figure 47: Intensity profiles for sinusoidal Laguerre-Gauss modes u^f . The Upo modes are identical 
to the helical modes. However, for azimuthal mode indices / > the pattern shows I dark radial 
lines in addition to the p dark concentric rings. 



to be computed (or at which coupling equations are to be defined). In our experience the quality 
of a computation or simulation and the correctness of the results depend critically on the choice 
of these beam parameters. One might argue that the choice of a basis should not alter the result. 
This is correct, but there is a practical limitation: the number of modes having non-negligible 
power might become very large if the beam parameters are not optimised, so that in practice a 
good set of beam parameters is usually required. 

In general, the Gaussian beam parameter of a mode is changed at every optical surface in a well- 



defined way (see Section 7.11 1. Thus, a possible method of finding reasonable beam parameters for 
every location in the interferometer is to first set only some specific beam parameters at selected 
locations and then to derive the remaining beam parameters from these initial ones: usually it is 
sensible to assume that the beam at the laser source can be properly described by the (hopefully 
known) beam parameter of the laser's output mode. In addition, in most stable cavities the 
light fields should be described by using the respective cavity eigenmodes. Then, the remaining 
beam parameters can be computed by tracing the beam through the optical system. 'Trace' in 
this context means that a beam starting at a location with an already-known beam parameter is 
propagated mathematically through the optical system. At every optical element along the path 
the beam parameter is transformed according to the ABCD matrix of the element (see below). 

7.11 ABCD matrices 

The transformation of the beam parameter can be performed by the ABCD matrix-formalism |361 
[5^ . When a beam passes an optical element or freely propagates though space, the initial beam 
parameter qi is transformed into (72 • This transformation can be described by four real coefficients 
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as follows: 



with the coefficient matrix 



92^ _ 

M = 



rii 



B 



C 91 



D' 



A B 
C D 



(154) 



(155) 



ni being the index of refraction at the beam segment defined by qi , and 712 the index of refraction 
at the beam segment described by q2. ABCD matrices for some common optical components are 
given below, for the sagittal and tangential plane. 

Transmission through a mirror: A mirror in this context is a single, partly-reflecting surface 
with an angle of incidence of 90°. The transmission is described by 



M 




Qi 



Q2 



Figure 48: 

with Rq being the radius of curvature of the spherical surface. The sign of the radius is defined 
such that i?c is negative if the centre of the sphere is located in the direction of propagation. The 



curvature shown above (in Figure 48), for example, is described by a positive radius. The matrix 
for the transmission in the opposite direction of propagation is identical. 

Reflection at a mirror: The matrix for reflection is given by The reflection at the back surface 



M 



1 

2ni 

Rc 




1 



92 



Figure 49: 

can be described by the same type of matrix by setting C = 2n2/i?c- 

Transmission through a beam spHtter: A beam splitter is understood as a single surface 
with an arbitrary angle of incidence ai. The matrices for transmission and reflection are different 
for the sagittal and tangential planes (Afg and Mt): with a2 given by Snell's law: 



and An by 



An 



ni sin (ai) = 712 sin (0:2), 

712 cos (a2) — ni cos (ai) 
cos (ofi) cos (q!2) 



(156) 
(157) 
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cos OL2 
COS Ctl 

An 

Rc 



1 

An 

-Rc 



cos a.1 

cos CK2 





^^2 























Figure 50: 



If the direction of propagation is reversed, the matrix for the sagittal plane is identical and the 
matrix for the tangential plane can be obtained by changing the coefficients A and D as follows: 



A 
D 



1/D. 



(158) 



Reflection at a beam splitter: The reflection at the front surface of a beam splitter is given 
by: To describe a reflection at the back surface the matrices have to be changed as follows: 



1 

2ni 



-Rc cos cKi 



1 

2ni cos ai i 

Rc ^ 



ni \ 


"2 























Figure 51: 



-Rc — > —Rc, 

ni — > 712, (159) 
ai — > —a2- 

Transmission through a thin lens: A thin lens transforms the beam parameter as follows: 
where / is the focal length. The matrix for the opposite direction of propagation is identical. Here 



Figure 52: 

it is assumed that the thin lens is surrounded by 'spaces' with index of refraction n — 1. 

Transmission through a free space: As mentioned above, the beam in free space can be 
described by one base parameter qq. In some cases it is convenient to use a matrix similar to that 
used for the other components to describe the z-dependency of q{z) = qq + z. On propagation 



A 



V 



(12 
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through a free space of the length L and index of refraction n, the beam parameter is transformed 
as follows. The matrix for the opposite direction of propagation is identical. 



M 



1 ^ 

n 

1 



n 



q2 



Figure 53: 
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8 Interferometer Matrix with Hermite— Gauss Modes 

In the plane- wave analysis Section |1.4[ a laser beam is described in general by the sum of various 
frequency components of its electric field 

i?(t,z)=^aj exp (^i{LUj t ~ kjz)^ . (160) 
j 

Here we include the geometric shape of the beam by describing each frequency component as a 
sum of Hermite-Gauss modes: 

E{t,x,y, z) = ^ Qj y) exp (i (wj t - fcjz)). (161) 

j n , m 

The shape of such a beam does not change along the z-axis (in the paraxial approximation) . More 
precisely, the spot size and the position of the maximum intensity with respect to the z-axis may 
change, but the relative intensity distribution across the beam does not change its shape. Each 
part of the sum may be treated as an independent field that can be described using the equation 
for plane- waves with only two exceptions: 

• the propagation through free space has to include the Gouy phase shift, and 

• upon reflection or transmission at a mirror or beam splitter the different Hermite-Gauss 
modes may be coupled (see below). 

The Gouy phase shift can be included in the simulation in several ways. For example, for reasons 
of flexibility the Gouy phase has been included in Finesse as a phase shift of the component space. 

8.1 Coupling of Hermite-Gauss modes 

Let us consider two different cavities with different sets of eigenmodes. The first set is charac- 
terised by the beam parameter qi and the second by the parameter q2. A beam with all power 
in the fundamental mode uoo(9i) leaves the first cavity and is injected into the second. Here, two 
'misconfigurations' are possible: 

• if the optical axes of the beam and the second cavity do not overlap perfectly, the setup is 
called misaligned, 

• if the beam size or shape at the second cavity does not match the beam shape and size of the 
(resonant) fundamental eigenmode ((?i(^cav) 7^ Q2{zciiv)), the beam is then not mode-matched 
to the second cavity, i.e., there is a mode mismatch. 

The above misconfigurations can be used in the context of simple beam segments. We consider 
the case in which the beam parameter for the input light is specified. Ideally, the ABCD matrices 
then allow one to trace a beam through the optical system by computing the proper beam parameter 
for each beam segment. In this case, the basis system of Hermite-Gauss modes is transformed in 
the same way as the beam, so that the modes are not coupled. 

For example, an input beam described by the beam parameter qi is passed through several 
optical components, and at each component the beam parameter is transformed according to the 
respective ABCD matrix. Thus, the electric field in each beam segment is described by Hermite- 
Gauss modes based on different beam parameters, but the relative power between the Hermite- 
Gauss modes with different mode numbers remains constant, i.e., a beam in a uqo mode is described 
as a pure uqo mode throughout the entire system. 
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In practice, it is usually impossible to compute proper beam parameter for each beam segment 
as suggested above, especially when the beam passes a certain segment more than once. A simple 
case that illustrates this point is reflection at a spherical mirror. Let the input beam be described 



by qi . From Figure 49 we know that the proper beam parameter of the reflected beam is 



with Rq being the radius of curvature of the mirror. In general, we get qi ^ q2 and thus two 
different 'proper' beam parameters for the same beam segment. Only one special radius of curvature 
would result in matched beam parameters [qi =92)- 



8.2 Coupling coefRcients for Hermite— Gauss modes 

Hermite-Gauss modes are coupled whenever a beam is not matched to a cavity or to a beam 
segment or if the beam and the segment are misaligned. This coupling is sometimes referred to as 
'scattering into higher-order modes' because in most cases the laser beam is a considered as a pure 
TEMqo mode and any mode coupling would transfer power from the fundamental into higher-order 
modes. However, in general, every mode with non-zero power will transfer energy into other modes 
whenever mismatch or misalignment occur, and this effect also includes the transfer from higher 
orders into a low order. 

To compute the amount of coupling the beam must be projected into the base system of the 
cavity or beam segment it is being injected into. This is always possible, provided that the paraxial 
approximation holds, because each set of Hermite-Gauss modes, defined by the beam parameter 
at a position z, forms a complete set. Such a change of the basis system results in a different 
distribution of light power in the new Hermite-Gauss modes and can be expressed by coupling 
coefficients that yield the change in the light amplitude and phase with respect to mode number. 

Let us assume that a beam described by the beam parameter qi is injected into a segment 
described by the parameter 52- Let the optical axis of the beam be misaligned: the coordinate 
system of the beam is given by (a;, y, z) and the beam travels along the z-axis. The beam segment is 
parallel to the z'-axis and the coordinate system (a;', y', z') is given by rotating the (x, y, z) system 
around the j/-axis by the misalignment angle 7. The coupling coefficients are defined as 

{qi) exp {i {ujt - kz)j = ^ ((J2) exp (^i (wt - fcz')^ , (163) 

n' ,ni' 

where Unm{qi) are the Hermite-Gauss modes used to describe the injected beam and Un'm'iq2) 
are the 'new' modes that are used to describe the light in the beam segment. Note that including 
the plane wave phase propagation within the definition of coupling coefficients is very important 
because it results in coupling coefficients that are independent of the position on the optical axis 
for which the coupling coefficients are computed. 

Using the fact that the Hermite-Gauss modes Unm are orthonormal, we can compute the 
coupling coefficients by the convolution |^ 

Since the Hermite-Gauss modes can be separated with respect to x and y, the coupling coefficients 
can also be split into k„mn'm' = knn'kmm'- These equations are very useful in the paraxial 
approximation as the coupling coefficients decrease with large mode numbers. In order to be 
described as paraxial, the angle 7 must not be larger than the diffraction angle. In addition, to 
obtain correct results with a finite number of modes the beam parameters qi and q2 must not differ 
too much. 
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The convolution given in Equation (164) can be computed directly using numerical integration. 
However, this is computationally very expensive. The following is based on the work of Bayer- 
Helms "B . Another very good description of coupling coefficients and their derivation can be found 
in the work of Vinet |55j . In [S] the above projection integral is partly solved and the coupling 
coefficients are given by simple sums as functions of 7 and the mode mismatch parameter which 
are defined by 

K =]^{K^ + iK2), (165) 

where Kq = (z/j — z'^/z'j^ and K2 = {{z — zq) — (z' — z'q))/z'^. This can also be written using 
g = i zr + z — zo , as 

The coupling coefficients for misalignment and mismatch (but no lateral displacement) can then 
be written as 

fc„„, = (-l)"'£;(")(n!n'!)i/2(l + i^o)"/2+i/^(l + if*)-("+"'+i)/2 {Sg - Su} , (167) 

where 

[n/2l In' /2] _ „ / o / min(u,u') _ / 

Sg^E E (n-2M)'(»'-2p')! ^ a (2a) ! (^ " a) ! (/^ - a) ! , 

[(n-l)/2l \(n'-l)/2] _ „ , / „ / , min(u,u') _ , 

— (n-2/j-l)!(n'-2Ai'-l)! ^ (2ct+1)!(^-o-)!(ai'-(t)! " 

fi=0 p'=0 a=0 



(168) 



The corresponding formula for fcmm' can be obtained by replacing the following parameters: n 
m, n! m! , X,X {) and E^^^ — > 1 (see below). The notation [n/2] means 



'm' 









ni/2 if m is even, 

(m — l)/2 if m is odd. 



(169) 



The other abbreviations used in the above definition are 



X^iiz'^ - z') sin (7)/(\/mFu;o), 
X = (izR + z')sin (7)/(Vl + if*u'o), 

F^K/{2{l + K^)), (170) 
F = K*/2, 

= cxp f-^^' 



In general, the Gaussian beam parameter might be different for the sagittal and tangential 
planes and a misalignment can be given for both possible axes (around the y-axis and around the 
a;-axis), in this case the coupling coefficients are given by 

where fc„„' is given above with 

and (172) 
wq wtfi, etc. 



71 



and 7 — > 7y is a rotation about the j/-axis. The k^m' can be obtained with the same formula, with 
the following substitutions: 



n — > m, 
n' m! , 

thus 

Wq Wsn, etc. 



(173) 



and 7 -T> 7a; is a rotation about the x-axis. 

At each component a matrix of coupling coefficients has to be computed for transmission and 



reflection; see Figure 54 



mm. 



Mm 



mm/A 

I'//////////////////. 



Figure 54: Coupling coefficients for Hermite-Gauss modes: for each optical element and each 
direction of propagation complex coefficients fc for transmission and reflection have to be computed. 
In this figure fci, fc2, fca, ^4 each represent a matrix of coefficients knmn'rn' describing the coupling 

of Un,m into Un' .rn' ■ 
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8.3 Finesse examples 
8.3.1 Beam parameter 




Figure 55: Finesse example: Beam parameter 

This example illustrates a possible use of the beam parameter detector 'bp': the beam radius of the 
laser beam is plotted as a function of distance to the laser. For this simulation, the interferometer 
matrix does not need to be solved, 'bp' merely returns the results from the beam tracing algorithm 
of Finesse. 

Finesse input file for 'Beam parameter' 

laser il 1 nl 7. laser with P=1W 

gauss gl il nl Im -2 '/, a dummy beam parameter 

maxtem 7, we need only the u_00 mode 

s sl 1 nl n2 7o a space of Im length 

bp width X w n2 7o detecting the beam width (horizontal) 

xaxis sl L lin 0.1 8 200 7. tuning the length of sl 



8.3.2 Mode cleaner 

This example uses the 'tem' command to create a laser beam which is a sum of equal parts in uqq 
and uio modes. This beam is passed through a triangular cavity, which acts as a mode cleaner. 
Being resonant for the uqq, the cavity transmits this mode and reflects the uio mode as can be 
seen in the resulting plots. 

Finesse input file for 'Mode cleaner' 

laser il 1 nl 7. laser with P=1W 

maxtem 1 7o need Hermite-Gauss modes up to n+m=l 

tem il 1 7« laser beam is a mix of u_00 and u_10 

tem il 1 1 

s sl 1 nl n2 7. a space of Im length 

7, triangular mode cleaner cavity 

bs bsl .9 . 1 n2 nrefl n3 n4 '/, input mirror 
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fexample_mc 
3i 



transmitted 



reflected 
3i 



Sat Aug 22 00:34:54 2009 



o 
o 



S, 01 



o 



2 
X [x/wxO] (ccd) 



o 
o 



S, 
^ I 



00 







X [x/wxO] (ccd) 
Figure 56: Finesse example: Mode cleaner 



s scl 2 n3 n5 '/, distance between bl and bs2 

bs bs2 .9 .10 ntrans dump n5 n6 7, output mirror 
s sc2 49 n4 n7 



s sc3 49 n6 n8 

bs bs3 1 n7 n8 dump dump 
attr bs3 Rc 150 
cav cavl bsl n3 bsl n4 
runl : beam ccd ntrans 
run2: beam ccd nrefl 
xaxis ccd x lin -3 3 200 
x2axis ccd y lin -3 3 200 
yaxis abs 



'/, distance between bl and bs3 
'/, distance between b2 and bs3 
end mirror 
Rc=150m for bs3 
'/o computing cavity parameters 
y, beam shape in transmission 
beam shape in reflection 
tuning x,y axes of beam detector 



'L plotting the absolute intensity 



8.3.3 LG33 mode 

Finesse uses the Hermite-Gauss modes as a base system for describing the spatial properties 
of laser beams. However, Laguerre-Gauss modes can be created using the coefficients given in 



Equation ( 149 ) . This example demonstrates this and the use of a 'beam' detector to plot amplitude 



and phase of a beam cross section. 
Finesse input file for 'LG33 mode' 

laser il 1 nl laser with P=1W 

gauss gl il nl Im a dummy beam parameter 

maxtem 9 we need modes up to n+m=9 

tem il HG coefficients to create LG33 mode 
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tem 
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164063 
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s si 1 


nl n2 







beam ccd n2 */, beam detector for carrier light 

xaxis ccd x lin -5 5 200 tune x position of beam detector 

x2axis ccd y lin -5 5 200 7. tune y position of beam detector 

yajxis abs:deg 7, plot amplitude and phase 

multi 
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A The Interferometer Simulation Finesse 



Throughout this document we have provided a number of text files that can be used as input files 
for the interferometer simulation Finesse [201 [23] • Finesse is a numerical simulation written in 
the C language; it is available free of charge for Linux, Windows and Macintosh computers and 
can be obtained online: http: //www.gwoptics . org/f inesse/' 

Finesse provides a fast and versatile tool that has proven to be very useful during the design 
and commissioning of interferometric gravitational- wave detectors. However, the program has been 
designed to allow the analysis of arbitrary, user-defined optical setups. In addition, it is easy to 
install and use. Therefore Finesse is well suited to study basic optical properties, such as, the 
power enhancement in a resonating cavity and modulation-demodulation methods. 

We encourage the reader to obtain Finesse and to learn its basic usage by running the included 
example files (and by making use of its extensive manual). The Finesse input files provided in this 
article are in most cases very simple and illustrate single concepts in interferometry. We believe 
that even a Finesse novice should be able to use them as starting points to play and explore 
freely, for example by changing parameters, or by adding further optical components. This type 
of 'numerical experimentation' can provide insights similar to real experiments, supplementing the 
understanding through a mathematical analysis with experience and intuitions. 
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